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1.0 INTRODUCTION 


The purpose of this document is to describe the Airport-noise 
jLevels and Annoyance Model (ALAMO) in terms of a detailed description 
of the constituent modules, to describe the execution of ALAMO procedure 
files, necessary for system execution, and to present the source code 
documentation associated with code development at Langley Research 
Center. The modules constituting ALAMO are presented both in flow graph 
form, and through a description of the subroutines and functions that 
comprise them. 
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2.0 SYSTEM DESCRIPTION 


The Airport-noise Levels and Annoyance Model (ALAMO) is a system of 
five processor modules which when executed in sequence evaluate the 
effect of aircraft noise on a surrounding airport community. The 
five processors contained in the baseline ALAMO system are: the 

Integrated Noise Model, the Contour Processor, the SITE II Phase 1 
Processor, the SITE II Phase 2 Processor and the Report Generator. 
These processors use, in addition to the INM and SITE II databases, 
ten user named files to produce noise exposure and impact reports, 
as well as noise impacted population distribution plots overlaid on 
an octant line compass rose. These files are named as follows in 
the ALAMO submit file: XXX is the Airport Auxiliary file and is 
input to the Contour Processor and the Report Generator; XXX001 is 
Airport Description Data and is input to the Integrated Noise Model; 
XXX002 contains airport statistics generated from the Integrated 
Noise Model; XXX003 is noise contour data generated from the 
Integrated Noise Model to serve as input to the Contour Processor; 
XXX004 is the SITE II Phase 1 input file generated from the Contour 
Processor; XXX005 is the plot vector file generated in the Contour 
Processor; XXX006 is used as a temporary accumulator for the output 
generated from both the SITE II Phase 1 and Phase 2 processors; 
XXX007 contains the complete output from both Phase 1 and Phase 2, 
and serves as input to the Report Generator; XXX008 contains the 
Noise Impact Summary Report tabulated by the Report Generator; XXXR 
contains the Noise Impact Summary Report associated with a 
countermeasure case executed subsequent to the baseline case. 
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Figure 1. -Baseline System Job Flow 
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The logical unit numbers assigned by a particular processor is 
illustrated in figure 1. The integrated Noise Model (INM) Processor 
reads as input the INM database and Airport Description Data. The 
INM database contains standard information on aircraft noise and 
performance for commercial, general aviation, and military aircraft. 
The Airport Description Data (named file XXX001) contains runway, 
track, approach profile and traffic mix data for any particular 
airport. The INM Processor reads these files, generates an airport 
statistics file (named file XXX002), and computes noise contour data 
for noise levels 65 to 85 decibels. The noise contour data file 
(named file XXX003) is then input to the Contour Processor. 

The Contour Processor reads, in addition to the contour data file, 
the Airport Auxiliary File (named file XXX) . This file contains the 
latitude, longitude, displacement and translation values, and a user 
designated population growth rate for a given airport. The Contour 
Processor then generates noise footprint contour plots and writes 
them to the plot vector file (named file XXX005) , and computes a 
sequential file (named file XXX004) containing noise contour 
coordinates which is input to the SITE II Phase 1 Processor. 
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The SITE II Phase 1 Processor reads, in addition to noise contour 
coordinates the SITE II database. The SITE II database is 
contained on magnetic tape and can be accessed from the volume 
serial numbers ND0775, ND1062, NC0773 and NC0938, respectively. The 
census data for the states Alabama through Mississippi (in 
alphabetical order) are contained on tape number ND0775; the census 
data for the states Missouri through Wyoming are contained on tape 
number ND1062; the block group data which describes a physical area 
surrounding a landmark is contained on tape number NC0773; and ■ the 
census tract data is contained on tape number NC0938. The zipcode 
data contained oh tape number NC0753 may be optionally processed. 
The data is retrieved (see fig. 8) before SITE II processing and 
written to the direct access file NSF1, NSF2 and NSF3 to be input to 
the Phase 1 processor via the logical files TAPE21, TAPE22 and 
TAPE23. Tape number NC0773 or NC0938 is selected for retrieval 
depending on whether the user wants to process block groups or 
census tract data. The census data for each state is retrieved from 
tape numbers ND0775 and ND1062 and copied to the direct access files 
MNFl, MNF2 and MNF3 to be input to the Phase 2 processor via the 
logical files TAPE31, TAPE32 and TAPE33. The census information 
retrieved is demographic and housing data obtained from the 1980 
census, updated to 1981; additional census data was obtained from 
the 1970 census. The Phase 1 Processor selects records containing 
population values that are within the noise level boundary computed 
in the Contour Processor. 
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The Phase 2 processor compiles the contour and demographic 
data, arranges it and writes it to a file (named file XXX006) . 
The contents of the file XXX006 are transferred to the named file 
XXX007; this final file is input to the Report Generator. 

The Report Generator reads the data file written in SITE II and the 
Airport Auxiliary File (named file XXX) . Ten reports are generated, 
the first eight reports contain noise exposure correlations for each 
octant in the compass rose; the ninth report contains demographic 
noise exposure values and the tenth report (written on the named 
file XXX008) is the Noise Impact Summary Report. 
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3.0 MODULE DESCRIPTION 


The baseline ALAMO system was assembled by linking four major 
subsystems together to comprehensively assess the noise impact on 
various airport communities. The system is a combination of 
software developed by outside interest, and software developed here 
at Langley Research Center (LaRC) to augment ALAMO for use on the 
Network Operation System (NOS). The Integrated Noise Model (INM), 
which computes noise exposure values, was developed by the Federal 
Aviation Administration; the Contour Plot Processor, which plots the 
contour noise footprints and generates contour descriptions in SITE 
II format, was developed at LaRC by Computer Sciences Corporation 
(CSC) ; the SITE II demographic program, which provides access to the 
United States Census database, was developed by CACI, Inc.; and the 
Report Program, which generates noise impact statistics was 
developed at LaRC by CSC. The following four sections describe 
the logic flow of each subsystem program illustrated by a tree 
diagram of source modules. 
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3.1 INTEGRATED NOISE MODEL 


The Integrated Noise Model program VERSUN6 is logically divided into 
four distinct areas of code as illustrated in figure 2, input, 
processing, output and general utility routines. The Airport 
Description Data Deck and the INM database are the input files to 
VERSUN6. The Airport Description Data Deck is read from logical 
unit 7 in subroutine READIN. Since Fortran input/output routines 
add the characters TAPE as a prefix to the logical unit to form the 
file name, TAPE7 is the file name assigned to logical unit 7. From 
this point on all of the logical units will be referred to by their 
file names. 


Input 

Refer to the input area of routines in figure 2. Runway and track 
data are read in subroutine RWVRD, track data is input beginning at 
entry point TRAKRD in RWYRD. Subroutine HELG is referenced by RWYRD 
to compute circular segment variables. Aircraft approach 
parameters, altitude data, wind information and noise parameters are 
read in subroutines APPTRD, ALTRRD, WINDRD and NOISRD respectively. 
Aircraft performance profiles are read in and stored in subroutine 
PROFD. The aircraft takeoff performance parameters are corrected by 
airport altitude and ambient temperature in subroutine TPROF. User 
defined replacement values for the tolerances used for checking 
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Figure 2 .-Integrated Noise Model Flow 
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contour looping conditions are read in subroutine TOLRD. The 
aircraft mix data describing the type, magnitude and arrangement of 
aircraft activity at an airport is read in subroutine MIXRD. Refer 
to the Int egra ted Noise Model Version 2 User's G uid e (Section_2) for 
a detailed description of the Airport Description Data Deck. The 
INM database contains aircraft noise and performance data for each 
aircraft type. If the user decides to enter new aircraft 
definitions a call is made to subroutine LOAD to read the INM data- 
base. In the event that this occurs, a call is made to subroutine 
INFORM to inform anyone reading the output that the user has 
modified the aircraft data in the database or has created his/her 
own aircraft definitions. 

A call is made to subroutine METFIX to initialize the noise metric 
variables. The metric identification (NEF, LDN, CNEL, LEQ, ASDS or 
DOSE) is read from the second record in The Airport Description Data 
Deck. A call is also made to subroutine CHEKTHR to compare the 
thrust noise settings in the database. This call is made if the 
user has opted to modify the database. 
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3.1.2 Processing 


Each noise level contour is processed point by point. Initially 
subroutine FIRST is called from VERSUN6 to take the first guess at 
the new contour point position. A call is made to subroutine GRADIE 
which computes the new point by determining the gradient. If the 
first point has been computed subroutine NEWPNT is called to compute 
the next contour point. If the metrics NEF, LDN, CNEL or LEQ are 
used then NEWPNT is referenced, if the metrics ASDS or DOSE are used 
then a call is made to subroutine NEWASDS to compute the next 
contour point. After the position of each point is computed a call 
is made to subroutine EXPOSE to calculate the noise exposure value 
at that point. The first call in EXPOSE is made to subroutine SIFT 
to determine which flights create a significant contribution to the 
overall noise exposure level. Noise exposure correction values for 
track geometry are computed in subroutine HBT. The shortest 
distance from a contour point to a ground track segment are computed 
to determine the correction values. The distance to curved and 
straight segments are computed in subroutines CURVE and STRAIT 
respectively. The corrections are then computed for curved segments 
in the function CONTF and for straight segments in the function 
CONTG. For comparison purposes a call is made to subroutine ACNOIS 
to compute the noise exposure level for one single flight from any 
combination of metrics. To determine the noise exposure levels 
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for the ASDS and DOSE metrics calls are made to subroutines ASD2 and 
ASD3 respectively. For the ASDS metric a call is made to the time 
history routine TIHISI to determine the times that the selected 
threshold noise levels are exceeded by one flight. Subroutine 
OVERLAY is referenced to overlay the flight performance profile over 
the ground track segments. The aircraft position coordinates, the 
ground attenuation and the noise level at a point is computed in 
subroutines PSCOO, EGA and AL respectively. 


3.1.3 Output 


The output area of the code writes two files for output, the Airport 
Statistics file on TAPE6 and the Noise Contour Data file on TAPE2. 
Subroutine HEADER prints the column headers for the contour 
information. The formats for PLINE, which is an entry point in 
HEADER, are also set up. The actual contour points and values are 
printed a line at a time at the entry point PLINE. The Noise 
Contour Data file is written to TAPE2 (named file XXX003) in 
subroutine TPGRM. Subroutine INSIDE is called to determine if point 
is within the defined contour or not, if not, that point is not 
written to the output file. Subroutine SKFIL is called to skip pass 
the end of file mark to start a new file. The entry points SKPPFF 
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skip forward one file or skip backwards one file. 


3.1.4 Utilities 


The system utilities are subroutines and function subprograms that 
are used to perform routine vector, sort and arithmetic operations. 
Vector operations such as addition, subtraction, dot product, vector 
to vector transfer, unit vector computations, vector magnitude 
computations and vector by scalar multiplication are performed by 
the functions VADD, VSUB, VDOT, VTRN, VUNT, VMAG and VSCL 
respectively. The subroutine ZERO is used to clear vectors. The 
sort utility S0RT2 sorts three arrays of numerical data in ascending 
order. The entry point S0RT3 in S0RT2 sorts two arrays in ascending 
order. Ihe subroutine S0RT4 sorts numerical data from four arrays 
or tables in descending order. S0RT4 is referenced from subroutine 
MIXRD to assemble the flight identification table and the associated 
equivalent operations for the noise metrics NEF, LDN, CNEL or LEQ. 
The degrees to radians conversion is performed by subroutine DGTRD; 
value normalization and sign determination are performed by the 
routines SGNR and SGNI respectively. 
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3.2 CONTOUR PROCESSOR 


The contour processor program AUPHS1 converts the Noise Contour Data 
file (named file XXX003) generated in the INM processor, from x-y 
coordinate points to a set of polygons for each noise level. These 
polygons are in a format compatiable for input to the SITE II Phase 
1 processor (see fig. 1). The subroutine INIT is called to read the 
Airport Auxiliary File (named file XXX) from TAPE2. This file 
contains the latitude, longitude, and optional north-south and 
east-west distances from the airport latitude and longitude. Axis 
translation data, axis length information and the airport title are 
also contained in the Airport Auxiliary File. Subroutine SYSTEM is 
called to issue error messages if the file doesn't process properly. 
The Noise Contour Data file on TAPE5 is read and evaluated to 
determine the number of records between the metric card and the x-y 
coordinate points. The header data and the x-y points for the noise 
contour are read from TAPE5 in subroutine RDCONTR and the axis 
translation is performed in subroutine TRNSLT (see fig. 3). The 
point reduction software is driven by subroutine SITE2. There is 
150 point limit for each contour. Subroutine NEWCRV is called to 
delete all unnecessary points below a given maximum on a noise 
contour. The area of the new contour is then computed to make sure 
that there is no more than a one percent difference in area between 
the old contour and the new reduced point contour. The plot vector 
file is then initialized in subroutine PLOT if this is the first 
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Figure 3. -Contour Processor Flow 
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contour processed. The entry point GETCON in the COMPASS routine 
JOBCON is referenced to sample job control register three. This is 
done to determine if this is a baseline case run. Plot information 
for the compass rose, plot boundaries and labels are written to the 
plot vector file in subroutine ROSE. 

The SITE II input file contains four types of card images containing 
either the characters 'A', 'B', 'C', or 'D' . These characters 
indicate the following: 


•A' 

Run Identification 



'B' 

Site Identification 



*C' 

Case Identification 



'D' 

Case Identification Continuation 



The reader is referred to the SITE II User's Manual, section 
a detailed description of these card images. After the plot 

IV for 
vector 

file is 

generated SITE II input file is written to TAPE4. 

The 

run 


number is written on card type 'A' and the airport name along with 
the airport latitude, longitude, the north-south displacement value 
and the east-west displacement value. The subroutine OCTANT is 
called to coordinate the polygon definitions. Subroutine OCTTAB is 
called from OCTANT to build a table of start and end points for each 
line segment in each octant. Subroutine IOCTNT is called from 
OCTTAB to determine which octant an x-y point lies in. After the 
start/end table is built a call is made to subroutine ENRICH if two 
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are made to the functions INTABLE and IOCTNT to determine if the 
ordinal is in the start/end table and which octant the x-y point 
lies in, respectively. If two adjacent points that are not in the 
same octant are not in adjacent octants subroutine ADDPTS is called 
to add points. Subroutine POLY is then called from OCTANT to build 
polygons from the line segment table. Subroutine DSTNCE is called 
from POLY to find the distance from start/end points to the origin. 
Subroutine INTSCT is referenced from DSTNCE to compute the line 
equation for each segment, determine the octant of each point to 
find the octant line they straddle and to find the point of 
intersection on the octant line. Upon return to TOLY 'C' type cards 
are written to TAPE4 containing the metric name, a table of LDN 
metrics and the area of the polygon in square miles. Subroutine 
NTHSTH is then called from POLY to convert x-y points from feet to 
statue miles and to write 'D' type cards containing there converted 
points to TAPE4. A final call is made to subroutine PLTPOL to plot 
the polygons given the x-y points and a scale factor. Execution 
control is then returned to program PHASE1. A final call is then 
made to BLDARR to extract information from the scratch unit TAPE4 
and build the SITE II input file on TAPE1 (named file XXX004) . If 
any deck contains more than 150 points then that deck is subdivided 
in this routine. Since each 'A' type card indicates the beginning 
of a new SITE II input deck an end-of-file mark is written to 
separate decks and the run counter is bumped. 
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3.3 SITE II DEMOGRAPHIC MODULES 


There are two SITE II processor modules, the phase 1 program MAIN1 
(see fig. 4) and the phase 2 program MAIN2 (see fig. 5). The phase 
1 module reads the SITE II input file (named file XXX004) and 
generates an intermediate file of selected SITE II records for 
sorting. The phase 2 module reads the sorted file and generates the 
input file (named file XXX006) . The information on XXX006 is 
aggregated on the named file XXX007, which serves as input to the 
Report Generator Module. 


Phase 1 


The phase 1 program MAIN1 initially calls subroutine INITAL to set 
default values for data item counts, maximum number of polygon 
points per run, population estimates, logical unit assignments, 
record length and block sizes for file structure. A call is then 
made to subroutine INPUT to read the SITE II input file (named file 
XXX004) and insure that there is a valid sequence of card types. 
Calls are initially made from INPUT to the subroutines WTCTRL and 
CDCBLK, to open the selected record file on TAPE1 and to encode 
holler ith blanks to blanks configured in "A" format, respectively. 
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Figure 4. -SITE. II Phase 1 Processor flow 
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Figure 5. -SITE II Phase 2 Processor Flow 
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The SITE II input file is then read a record at atime and the first 
character on the record is checked for the card type. A call is 
then made to subroutine PROCES to define the structure for each card 
type. Calls are then made from PROCES to subroutines FREFOR and 
EDIT, to convert free form strings to fixed format card images and 
to edit each record to insure that each field has valid characters, 
respectively. Subroutine VERIFY is called to check specific data 
items in the SITE II input file. A call is then made to subroutine 
PUTDAT to create a sort key for the selected record file (TAPE1) . 
Execution control is then passed back to subroutine INPUT to read 
another input record. 

After the SITE II input has been read a call is made to subroutine 
GETZIP from MAIN1. This subroutine takes in array of requested zip 
codes and searches the zip file. For each zip found tract inclusion 
records are written to the selected record file on TAPE1. 
Subroutine RDZIP is called from GETZIP to open, read, rewind or 
close the zip file on TAPE1. A call is then made to subroutine 
LIMITS to find the latitude and longitude limits for each site and 
each run. Subroutine CREBND is referenced from LIMITS to convert 
sectors, chords and polygons to a form useable by subroutine SEARCH. 
Subroutine SEARCH reads the north/south records to determine the 
cases for which the record is a hit. For each hit a control record 
is written to TAPE1. Calls are made to subroutine NSPOSN, RDNTOS 
and NTCTRL, to get the first record below a given latitude, to get 
the next record and to write to the control file, respectively. 
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3.3.2 Phase 2 


Before the phase 2 processor is invoked, the track inclusion control 
file written by the phase 2 processor on TAPE1 is sorted using the 
CDC utility SORTMRG . The sorted control file is written on TAPE2 
for input to the phase 2 processor. Subroutine INITAL is referenced 
to initialize operational variables the same as the phase 1 
processor. Subroutine ACCUML is then called to accumulate the area 
data for each case. Subroutine RDCTRL is called from ACCUML to read 
the control file which contains area information for each site. 
Subroutine MNFIND is called to read the main file. The file header 
is read to check for the state and to determine' the number of 
records for that state. Upon return, ACCUML adds the percentage of 
area data to the total and writes the component area to TAPE3. 

There are three levels of data on the SITE control file; level 1 is 
RUN data which is the population summation of all the areas formed 
around the sites under study; level 2 is SITE data, this is the 
reference point around Which one or more boundaries are defined; and 
level 3 is CASE data which represents a physical area around the 
site. Subroutine OUTA controls the production of growth, profile 
and component area report data by referencing subroutine GETDAT to 
read the specified level of data from TAPE 2. A call is then made to 
subroutine CNTROL to resolve specific requests for output. Calls 



made to COMPON to generate values for all the reports, accumulate 
data for SITE and RUN totals, to write summary report data in 
condensed form, to produce profile report data and to write com- 
ponent area report data, respectively to TAPE6 for input to the 
Report Generator. Upon return to MAIN2 a final call is made to 
subroutine OUTB to control the production of columnar and comparison 
reports data through calls to subroutine COLUMN and CQMPAR. This 
data is also written to TAPE6 for input to the Report Generator. 
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3.4 REPORT GENERATOR 


The SITE II report file (named file XXX007) generated following the 
SITE II phase 2 processor is input to the Report Generator on TAPE5. 
The Airport Auxiliary File (named file XXX) is also input on TAPE1 

to the Report Generator program AUREPT (see fig. 6) to obtain the 

year of population, the population growth rate and the airport name 
and address. The SITE II case title is then read from the first 
record on the SITE II report file. The first call is made to 
subroutine BLDTAB to build a 3-dimensional table of sixteen 
demographic variables for each octant of each LDN band. The 

starting LDN index is defined by calling the entry point GETCON (in 

subroutine JOBCON) to fetch the job control register R3. A call is 
made to subroutine RDREPT to read the following demographic 
information from the SITE II report file: 

o Male and female age groups 
o Family income, home values and occupations 
o The number of single dwellings households 
o Households with T.V.'s and air-conditioners 

This information is stored in the labelled common block /COMM/ and a 
return is made back to subroutine BLDTAB. Beginning at the LDN value 
determined from job control register R3, the 3-dimensional table 
is generated, stored in the array TABLE and passed back to REPORTS 
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Figure 6. -Report Generator Flow 
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to print the heading for each octant. The TABLE array is checked for a 
void octant and a call is made to subroutine CHART. Subroutine CHART 
prints sixteen demographic variables by LDN band values and calls 
subroutine BREAK to transfer seven columns out of the 3-dimensional 
table to a vector array called LINE. An alternate call is made to 
subroutine CHART1 for special printing if an LDN contour is missing. A 
call is then made to subroutine SUMARY to print the Noise Impact Summary 
Report. The functions INTERP and PERCNT are referenced by SUMARY to 
compute the equivalent noise level value and the weighting factor for a 
given LEW value. Finally subroutine LASTPAG is referenced to write a 
separate Noise Impact Summary to TAPE2 (named file XXX008) . 
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■4.0 SYSTEM USAGE 

The baseline ALAMO system can be used by submitting two submit 
files, the SITE II file retrival file SITEJOB and the baseline ALAMO 
execution submit file BASJOB. These submit files along with the 
procedure files they invoke will be described in the following sections 
to give an overview of the logic flow of the system. 
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4.1 SITE II DATABASE RETRIEVAL 

The SITE II database resides on five magnetic tapes. These tapes 
are organized in the following manner: 


Volume Serial No. 

Tape Description 

ND0775 

State Code less than or equal to 28 

ND1062 

State Code greater than 28 

NC0773 

Physical area surrounding a landmark 

NC0938 

Census Tract Data 

NC0753 

Zip codes 


The submit file illustrated in figure 7 is used to retrieve the SITE II 
database. The procedure file SETUP which resides on the permanent file 
SET is invoked by BEGIN statement. The state integer codes 17, 21 and 
29, the state character codes IL, KY and MO along with the AREAS and 
ZIPS parameters are passed to the procedure file SETUP. In this example 
AREAS=1 indicates that the tape containing physical areas will be 
retrieved instead of the census tract data, ZIPS=0 indicates that the 
zip codes tape will not be retrieved. The procedure file SETUP is 
illustrated in figure 8. The parameters passed to SETUP through the 


PROC statement are as follows: 



A = 

17 

E = 

KY 

B = 

21 

F = 

MO 

C = 

29 

G = 

1 

D = 

IL 

H = 

0 
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DELIVERY Irro 


SJBFIL, TIES, 01100003. 

USER 

CHARGE .. ,LRC. 

GET, SET. 

KGIM, SETLP, SET, 17, 21 , 29, IL, KY, MO, 1 , 0. 

LaI 1 • 


EDITING NOTES: 

- REPLACE INTEGER CODES 17, 21 
APPROPRIATE INTEGER CODES FOR 
INTEREST 


AND 29 WITH ThE 
THE STATES OF 


- REPLACE CHARACTER CODES IL, 
APPROPRIATE CHARACTER CODES 
INTEREST 


KY, AND MO WITH TVE 
FOR TVE STATES OF 


Figure 7. -Submit File for SITE II Retrieval 
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. PROC, SOUP, A, B, C, D. E, F, G, H. 

BEGIN, T, SET, A, A, ND0775, ND1062, D, MSF1 . 
BEGIN, T, SET, fi,B,ND0775,ND1062,E,rTf r 2. 
BEGIN, T, SCT, B, C, ND0775, ND1062, F, rTF3 . 
RETURN, TAPE1 , fTF 1 , fTFE , fTF3 . 
IFE,(G.rE.0),LAEELA. 

VSN,TAPE2-NC0773. 

BEGIN, BT, SET, BLXGRP, D, NSE1 . 

BEGIN, BT , SET, BLKGRP, E, NSF2 . 

BEGIN, BT , SET , BLXGRP, F , NSF3 . 
SKIP(LABELB) 

ENDIF,LABELA. 

VSN,TAPE2-NC0938. 

BEGIN, BT, SET, TRACKS, D, NSF 1 . 

BEGIN, BT, SET, TRACKS, E, NSF2. 

BEGIN, BT, SET, TRACKS, F, NSF3. 
ENDIF,LABE1_B. 

RETURN, TAPE2. 

IFE, (H. EQ. 1 ) , LABELC. 

VSN, TAPES -NC0753. 

LABEL, TPPE3, PO-R, SI-2IPCDS, FI -ZIPALL. 
COPYBF, TRPE3, ZIPF 1 . 

REPLACE, Z I FF1. 

ENDIF, LABELC. • 

REVERT. MYCCL SETUP 


Figure 8. -SITE II Procedure File SETUP 



I 


1 


Since at least three states from the database must be retrieved for 
an ALAMO run, the procedure file T is invoked three times as indicated 
by the three BEGIN, T, ... statements. In the first statement the integer 
code 09 is passed twice through the A and B parameters and the state 
character code CT is passed through the D parameter. Since all of the 
database tapes are multi-file files containing a file per state, the D 
parameter is the file identifier (FI) parameters on the LABEL card. The 
procedure file T illustrated in figure 9 is invoked by the BEGIN 
statement and the A and B parameters are checked to determine which tape 
to load. If the state integer code is less than or equal to 28 then the 
tape with the volume serial number ND0775 is retrieved, if the integer 
code is greater than 28 then the tape numbered ND1062 is retrieved. The 
file containing data for the state of Connecticut (integer code 09, 
character code CT) is read from ND0775 and written to the direct access 
file MNFl; the file for the state of New Jersey (integer code 24, 
character code NJ) is read from ND0775 and. written to the direct access 
file MNF2; and the New York file (integer code 26, character code NY) is 
read from ND0775 and written to the direct access file MNF3. The G 
parameter is checked to determine if the block group or the census track 
tape will be read. In this example G = 1 indicates that the block group 
tape (NC0773) will be retrieved. Notice that the IFE statement does not 
work like an IF statement. If the expression is true, the statements 
following the IFE are executed until an ELSE or ENDIF with a matching 
label is encountered. The procedure file BT (see fig. 10) is invoked 
three times to write block group data from the states, Connecticut, New 
jersey and New York to the direct access files NSF1, NSF2 and NSF3. 


. PROC, T, A, B, T1 , T2, C, FT. 

PURGE, FF/THA, ST -LPF. 

DEFirC,FF. 

IFE,(A.LE.2B.AND.B.L£.28),LABELA. 
LABEL, TAPE1,PCHR,VSN-T1, SI-FNl, FI-C. 
SKIP, LABEL D. ~ 

ENDIF, LABELA. 

IFE,(A.GT.28.AND.B.GT.28),LAB£LB. 
LABEL, TAPEl,PO-R, VSN-T2, SI -rTS2, FI-C. 
SKIP,LAEELD. 

ENDIF, LABELS. 

RETURN, TAPE1. 

IFE,(B.GT.28),LABELC. 

LABEL, TAPE1 , PO -R, VSN-T2, SI - m2, FI -C . 
SKIP,LABELD. 

ENDIF, LABELC. 

LABEL, TAPE1 , PO-R, VSN-T1 , SI -mi , FI -C. 
ENDIF, LABELD, 

COPYBF, TAFE1,FF . 

RETURN, FF. 

REVERT. MYCCL T 


Figure 9. -SITE II Procedure File T 



.PR0C,BT,SID,X,MS1. 

PURGE, MSI /ISP, ST -LPT. 
DEFirC.MSl. 

LABEL, TRPE2, SI -SID, FI -X, PO-R. 
C0PYBF,TPP£Z,MS1. 

RETURN, MSI. 

REVERT. MYCCL BT 


Figure 10. -SITE II Procedure File BT 
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The reader should note that the direct, access files MNF1, MNF2, MNF3, 
NSF1, NSF2 and NSF3 must be created on the WFS cluster. The H parameter 
is checked to determine whether the zip codes tape (NC0753) is to be 
retrieved. If H = 1 then the tape file is written to the indirect 
access file ZIPl. After these files have been created a baseline ALAMO 
case, can now be executed. 
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4.2 ALAMO JOB t'LUW 

To- run a baseline ALAMO case the submit file BASJOB (see fig. 11) 
must be submitted. Due to the large amount of execution time required 
the user should submit BASJOB to a CYBER 175 (R or T) machine. The job 
control register R3 is set to zero indicating a baseline run. 

4.2.1 Integrated Noise Model 

The procedure file A1 is invoked using the BEGIN statement. There 
are 10 parameters naming the ALAMO system files. The XXX should be 
replaced with the F.A.A. airport identifier. The procedure file A1 
illustrated in figure 12 describes each subsequent procedure file 
invoked and the parameters used. The direct access file CONTORA 
containing the absolutes for the INM processor program VERSUN6 is 
attached. The INM data base DATABIN is also attached and copied to 
TAPE7 for processing. The Airport Description Data file (named file 
XXX001) is loaded from the permanent file A. Fiyure 13 illustrates the 
execution of the INM processor. The local file A is the Airport 
Description Data file input to the INM processor CONTORA on TAPE5, the 
parameter B is the Airport Statistics file (named file XXX002) output 
from CONTORA on TAPE6 and the Noise Contour Data file (named file 
XXX003) generated on TAPE2 is copied to the local file C for input to 
the Contour Processor. 
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DELIVERY IhFO 


BRSJ0e,T12SG3,Cn2733a3. 

USER, 

CHPRGE, ,LRC. 

VSfK TPFE51 -mYTYY, TAPE52- NZZZZZZ ) 

SETCR3-0) 

BEGIN, Pl,R2,XXX231,XXX£^,XXXKG,XXXE34,XXXt^,XXXE^,XXXB^,XXX333,XXXR >OC< 
REWIND, TRPE6. 

ROUTE, TfiFES, DC-LP. 

EXIT. 


EDITING NOTES: 

- INSERT fCTUPL REEL NUMBERS IN VSN STATETENT 

- REPLACE XXX WITH F.A.A. IDENTIFIER 


Figure 11. -Submit File for Baseline Execution 
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1 


PR0C,A1,A,B,C.D,E:,F,G,H,J,K. 

•y fggyyy 


c x cyrmr m a 


FILE A2 CCNTAIhE T>C 5 FOLLOWING PROCEDURE. FILES: 


fil 

Bl 

Cl 


D1 

G1 


DRIVER FOR Ft COTFEETE ALAMO RUM 
EXECUTES PHASE 2 OF ALAMO 
URITES FILES PRODDUCED BY ALAMO RUN TO 
TAPE ( TAPES 1 FOR SUCCESSFUL RUN, TAPE52 
FOR RUN WITH ERRORS) 

INITIALIZES TAPES USED IN ALAMO RUN 
GETERATES POPULATION GRID AND STANDARD 
TRACKS 


* 

* 

* 

* 


A1 PROCEDURE PARfirCTER IDENTIFICATION: 


A 

B 

C 

D 

E 

F 

G 

H 

J 

K 


INI IN=UT CASE 
LIST ABLE INI OUTPUT 

INI OUTPUT - CONTOUR POINTS (PHASE 2 If-PUT) 
PHASE 1 OUTPUT (SITE II IFFUT) 

SAVPLT FILE - OUTPUT FROM PHASE 1 
SITE II OUTPUT (PARTIAL) 

SITE II OUTPUT (COrPLETE) - INPUT FOR PHASE 2 
10TH PAGE ALAMO REPORT (BASELIhE CASE) 

10TH PACE ALAMO REPORT ( C(XJNTERMEASURE CASE) 
ALAMO AUXILIARY FILE (ILPUT FOR PHASES 1 AND 


2 ) 




* 

* IF R3-I, NOISE REDUCTION SIMULATION RUN IS IN PROGRESS; 

* A MODIFIED VERSION OF INI DATABASE ALREADY EXISTS 

* 

ATTACH, CONTORA. 

I FE, (R 3il€.l),LABELl. 

ATTACH, DATRBIN. 

COPYEI , DATA3IN, TAPE7 . 

REWIND, TAPE7. 

RETURN, DATRBIN. 

ENDIF,LABEL1. 

GET, A. 

* 

* EXECUTE THE INI MODULE OF ALAMO 

* 

CONTORA, A, B. 

SEND,B,DC-LP. 

RETURN, CONTORA. 

REWIND, TAFE2. 

COPYEI, TAPEH,C. 

SKIP(LABELA) 

EXIT. 

SKIP(LA5EL99) 

ENDIF(LABELA) 

RETURN, TAFE5, TAPES, TAPE2, TAFEB. 

RETURN, TAPES, TAPE3, TAPE7. 


Figure 12. -Baseline Procedure File A1 
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ATTACH, CCtiTCRA. 

I FE, (R 3.HE.1),LAEEL1. 

ATTACH, DATABIN. 

COPYEI , DATABIN, TAPE7 . 

REWIND, TAFE7. 

RETURN, DATABIN. 

ENDIF,LABEL1. 

GET, A. 

3k 

* EXECUTE THE Iltl MODULE OF PLPHO 

)k 

CONTORA, A,B. 

SEND,B,DC-L°. 

RETURN, CONTCRA. 

REWIND, TAPE2. 

COPYEI, TAPE2,C. 

SKIP(LABELA) 

EXIT. 

SKIP(LABEL99) 

ENDIF (LABELA) 

RETURN, TAPE5, TAPES, TAPE2, TAPEB. 

. RETURN, TAPE9,TAPE3,TAPE7. 

3k 

3k EXECUTE PHASE 1 

3k 

REWIND, C. 

GET, AUFHS1A,K. 

ALFVS1A,C, , , , ,K. 

HFP OFF. 

IFE* ( R3 ! EQ . 0 ), LABELS . 

R£FLAGE,K. 

SET CR1G-R2) 


Figure I3.-INM and Contour Processor Execution 
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4.2.2 Contour Processor 


The indirect access file AUPHS1A containing the absolutes for the 
Contour processor program PHASE1 is loaded along with the Airport 
Auxiliary File (named file XXX) . The noise contour data (named file 
XXX003) is input through the C parameter to the Contour processor 
AUPHS1A, along with the Airport Auxiliary File through the K parameter. 
A plot vector file SAVPLT is generated to plot the noise polygons. The 
SITE II (named file XXX004) input file is also generated on TAPE1 and 
copied to the D file to serve as input to the SITE II phase 1 processor. 

4.2.3 SITE II Demographic Processor 


Prior to running the SITE II phase 1 processor, the direct access 
files containing the block group or census track portion of the SITE II 
database are attached as logical unit files TAPE21, TAPE22, and TAPE23 . 
The indirect access file PH1A containing the SITE II phase 1 processor 
program MAIN1 is loaded. The SITE II input file is input to PH1A 
through the D parameter and the input reflection file (named file 
XXX006) is output through the F parameter after execution. The SITE 
II control file is generated on TAPE1 to be sorted using the CDC sort/ 
merge utility SORTMRG. 
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The direct access files containing the state census data portion of the 
SITE II database are attached as logical unit files TAPE31, TAPE32 and 
TAPE33. The indirect access file PH2A containing the SITE II phase 2 
processor program MAIN2 is loaded and executed. The SITE II report file 
(named file XXX007) is output through the F parameter. The p file is 
then copied to the G file for input to the Report Generator. 

4.2.3 Report Generator 


The procedure file B1 is invoked using the BEGIN statement to 
process the Report Generator (see fig. 15). The listing of B1 describes 
the parameters used in the procedure file as illustrated in figure 16. 
The indirect access file AUREPTA containing the Report Generator program 
REPORTS is loaded, along with the Airport Auxiliary File (named file 
XXX) . The SITE II report file is input to AUREPTA through the G 
parameter and the Airport Auxiliary File is input through the k 
parameter. After execution, the noise exposure correlation and 
demographic noise exposure reports are written on TAPE6 and the 10th 
page ALAMO report is written and saved on the permanent file j (named 
file XXXR) . The reports on TAPE6 are printed before job control is 
returned to the procedure Al. 
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PR0C,B1,G,H,J,K. 


* B1 EXECUTES PHASE 2 OF ALAMO 

* 

* B1 PROCEDURE PARAMETER IDENTIFICATION: 

* G : SITE II REPORTS (PHASE 2 ItFUT) 

* H : 10TH PAGE ALAMO REPORT (COUNTERMEASURE CASE) 

* j : ierm page alamo report (baseline case) 

* K : ALAMO AUXILIARY FILE 


GET, ALREPTfl. 

GET,K. 

REWIND, G. 

PACX,G. 

f-pp opy, 

PLIREPTA^G, TAPES, , ,K,H. 
REWIND, H, TAPES. 

COPYEI, TAPES. 

IFE, (R3.EQ.0),LAEEL4. 
RE^TRCE, G. 

COPYEI, H,J. 

REWIND, J. 

SAVE J. 

REWIND, H,K,TAPES. 
PACX,K. 

REPLACE,*. 

ELSE, LA3EL4. 

APPEND, J,H. 
ENDIF,LABEL4. 

REVERT. NYCCL 


Figure 14 .-Report Generator Execution 
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PR0C,B1,G,H,J,K. 


* B1 EXECUTE5 PHAS E 2 OF ALPflO 


* 


* B1 PROCEDURE PmRPTETER IDENTIFICATION: 

* G : SITE II REPORTS (PHASE 2 IhFUT) 

* H : 10TH PAGE PLATO REPORT ( COLTfTERrEPiSLRE CASE) 

* J : 1ETTH POSE PLPDO' REPORT (BASELINE CASE) 

* K : PLPDO AUXILIARY FILE 


GET,ALFEPTR. 

GET,K. 

REWIND, G. 

PACK,G. 

DPP, OFF. 

alrepta.g, tapes, , ,k,h. 

REWIND, H, TAPES. 

COPYEI , TAPES. 

IFE, (R3.EQ.0),LABEL4. 
REPLACE, G. 

COPYEI, H,J. 

REWIND, J. 

SAVE J. 

REWIND, H,K, TAPES. 
PACX,K. 

REPLACE, K. 

ELSE , LABEL 4 . 

PPPEND, J,H. 
ENDIF,LABEL4. 

REVERT. MYCCL 


Figure 15 .-Baseline Procedure File B1 
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APPENDIX A 


SOURCE DOCUMENTATION 

The following listings are the source code documentation for the 
auxiliary programs developed here at LaRC for the ALAMO System. The 
documentation consist of two program modules and their corresponding 
subroutines and function subprograms; the contour data processor program 
AUPHS1 and the report generator program AUREPT. Each routine 
description contains the objective, the method used (where applicable), 
the input/output parameters, the calling routines, the routines 
referenced, a local variable dictionary, a description of common block 
variables and a description of the files used. 
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PROGRAM AUPHS1 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OBJECTIVE : 

THIS PROGRAM PROCESSES THE INTEGRATED NOISE MODEL (INM) 
CONTOUR OUTPUT FILE. THE MAIN OBJECTIVE IS TO CONVERT 
THE X, Y CONTOUR DATA INTO A SET OF POLYGONS IN AN INPUT 
FORM WHICH IS COMPATIBLE WITH SITE II, A DEMOGRAPHICS 
PROGRAM. AS THE POLYGONS ARE DEFINED THEY ARE PLOTTED 
USING A STANDARD X,Y PEN-PLOTTER (CALCOMP) . THE PROGRAM 
ALSO PERFORMS SOME SECONDARY FUNCTIONS AS REQUIRED bY TOE 
LIMITATIONS OF THE PROGRAMS THAT ARE DEPENDENT ON THE 
DATA IT PRODUCES. A NOTABLE SECONDARY FUNCTION IS THE 
REDUCTION OF DATA POINTS ON THE INM CONTOURS TO CONFORM 
TO THE SITE II RESTRICTION OF NO MORE THAN 15U POINTS PER 
SITE II EXECUTION. 


METHOD: 

THE PROGRAM READS THE X,Y DATA REPRESENTING A CONSTANT 
LEVEL NOISE METRIC CONTOUR AND REDUCES THE NUMBER OF X,Y 
POINTS TO LESS THAN 150. AN OCTANT LINE COMPASS ROSE IS 
SUPERIMPOSED OVER THE CONTOURS AND POLYGONS ARE DEFINED 
FOR EACH OCTANT BY TOE INTERSECTION OF THE OCTANT LINES 
AND THE CONTOUR LINE. TOE POLYGON X,Y POINTS ARE 
CONVERTED TO AN 80 CHARACTER PER LINE INPUT DECK THAT 
MEETS THE INPUT SPECIFICATIONS FOR SITE II. THE X,Y 
POLYGONS ARE ALSO PLOTTED VIA TOE LRCGOSF GRAPHICS 
ROUTINES. SUPPORTING DATA SUCH AS AIRPORT CENTROID 
LATITUDE AND LONGITUDE AND AIRPORT NAME AND ADDRESS ARE 
READ IN FROM AN EXTERNAL FILE NAMED TAPE2 . 


SUBPROGRAMS USED: 
BLDARR : 

CALPLT : 

CHKREC : 
IN IT : 

OCTANT : 

PLOT : 


RDCONTR: 

RESTPLT: 


EVALUATES TOE SITE II CARDS FOR 150+ POINTS 
VIOLATION AND ADJUSTS IF NECESSARY. 

PLOTTER PEN CONTROL ROUTINE; USED HERE TO 
TERMINATE THE PLOT FILE. 

LOOKS tOR A MESSAGE ON TOE CONTOUR INPUT FILE. 
READS THE SUPPLEMENTAL DATA FILE, TAPE2 AND 
INITIALIZES CERTAIN CRITICAL VARIABLES. 

DRIVES THE CONSTRUCTION OF AN OCTANT ORDINAL 
TABLE AND TOE DEFINITION OF TOE POLYGONS. 
INITIALIZES THE PLOT VECTOR FILE, COMPUTES THE 
PLOT SCALE FACTOR AND CAUSES TOE COMPASS ROSE 
TO BE DRAWN. 

INPUTS TOE DATA FOR 1 NOISE CONTOUR. 

ROUTINE USED TO PLOT THE RANGE OF THE LDN VALUES 
FOUND 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

€ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SETCON : COMPASS ROUTINE ENTRY POINT FOR SETTING THE 3 

JOB CONTROL REGISTERS. 

SITE2 : DRIVES THE POINT REDUCTION PROCESS. 

TRNSLT : PERFORMS X,Y POINT TRANSLATION ON X,Y CONTOUR 

POINTS BASED ON DATA ON FILE TAPE2. 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

EWDIS : EAST/WEST DISPLACEMENT FROM LATITUDE POINT. 

FACT : THE LENGTH OF THE POSITIVE X AND Y AXIS IN 

INCHES. 

H : X DISTANCE TRANSLATION. 

ICASE : STORAGE FOR CASE NAME INFORMATION. 

ISW : INTEGER VARIABLE USED AS AN INDICATOR OF RECORD 

STRUCTURE OF TAPE5 

ITITLE : 8 WORD ARRAY STORING THE AIRPORT TITLE. 

JJ : INTEGER USED AS LOOP COUNTER 

K : Y DISTANCE TRANSLATION. 

KEOF : END-OF-FILE FLAG SET TO 1 BY RDCONTR TO 

INDICATE E-O-F ON CONTOUR INPUT FILE. 

LAT : LATITUDE OF AIRPORT CENTROID IN CHARACTER FORM. 

LDN : NOISE METRIC VALUE IN CHARACTER FORM. 

LINE : CHARACTER STORAGE FOR NOISE METRICS AS THEY 

ACCUMULATE. 

LONG : LONGITUDE OF AIRPORT CENTROID IN CHARACTER FORM 

MET : NOISE METRIC IN CHARACTER FORM; EG: "LDN" 

NCON : 18 BIT INTEGER VALUE TO SET JCR R1 VIA SETCON. 

NEWPTS : THE NUMBER OF X,Y POINTS AFTER POINT REDUCTION. 

NPTS : THE NUMBER OF X,Y POINTS BEFORE POINT REDUCTION 

NR3 : JOB CONTROL REGISTER, R3. 

NSDIS : NORTH/SOUTH DISPLACEMENT FROM LONGITUDE POINT. 

NSTCON : 18 BIT INTEGER VALUE TO SET JCR R2 VIA SETCON. 

RUNCNT : COUNTER OF CONTOURS CORRESPONDING TO ONE SITE 

RUN PER COUNT OF RUNCNT. 

SCALE : THE PLOTTING SCALE FACTOR. 

X : THE X POINTS OF THE CONTOUR BEFORE POINT REDUCT 

XNEW : THE X POINTS OF THE CONTOUR AFTER POINT REDUCT. 

Y : ■ THE Y POINTS BEFORE REDUCTION. 

YNEW : THE Y POINTS AFTER REDUCTION. 


COMMON BLOCK CROSS-REFERENCE LIST: 
NONE 
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c 

c 

c 

c 

c 
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c 

c 
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c 
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c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMMON VARIABLES: 

/ITITL/ 

ICASE1 : 4-WORD ARRAY CONTAINING CASE TITLE 

/PLTINFO/ 

CENT : THE CENTER X AND Y VALUE OF THE COMPASS 
ROSE 

HT : THE HEIGHT IN INCHES OF. CHARACTERS TO BE 
NOTATED 

ISTRIN : CHARACTER STRING (2 WORDS) USED TO NOTATE 
RANGE OF LEW VALUES 

S : THE DISTANCE FROM THE LEADING EDGE OF A 

SYMBOL TO THE LEADING EDGE OF THE NEXT 
SYMBOL 

W : REAL VARIABLE REPRESENTING THE WIDTH OF 

ONE SYMBOL 


FILES USED: 

INPUT: 

TAPE5 - SEQUENTIAL INPUT FILE CONTAINING ThE CONTOUR 
X,Y POINTS AND SUPPORTING DATA IN DISPLAY 
CODE FORMAT. 

CARD1 DATE IN YY/MM/DD. FORM. 

CARD2 INM BANNER AND VERSION ID. 

CARD3 NOISE METRIC AND VALUES. 

CARD4 CASE TITLE 

CARD5 NULL RECORD OR SPECIAL MESSAGE. IF A 

SPECIAL MESSAGE EXISTS, IT IS FOLLOWED 
BY A VARIABLE NUMBER OF EMPTY RECORDS 
(CARD5B, CARD5C, ETC.) . 

CARD6 6111, 6E11.5 WHERE THE 6TH INTEGER IS 
THE NUMBER OF CONTOUR POINTS. 

CARD7-16 (10E14.7 FORMAT) THE 1ST 10 X POINTS. 
CARD17-25 THE 2ND 10 X POINTS. 


CARD K-K+9 (10E14 .7 FORMAT) THE 1ST 10 Y POINTS. 


E-O-F 

NEXT CONTOUR AS ABOVE 


A- 4 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 
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c 

c 

c 
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TAPE2 - SEQUENTIAL INPUT FILE CONTAINING 

SUPPLEMENTAL DATA NOT GIVEN ON TAPE5. 
DISPLAY CODE FORMAT: 

CARDS 1-6 - INFORMATION USED BY OTHER ALAMO 

PROGRAMS THAT SHARE THIS FILE 
CARD7 LATITUDE, LONGITUDE, NORTH/SOUTH AND 

EAST/WEST DISPLACEMENT. 

CARD8 H, K TRANSLATION VALUES. 

CARD9 POSITIVE X,Y AXIS LENGTH tOR PLOT. 

CARD10 AIRPORT TITLE 

CARD11 AIRPORT ADDRESS. 

OUTPUT: 

TAPE6 - CONVENTIONAL FORTRAN STANDARD OUTPUT FILE; 
CURRENTLY NOT USED. 

TAPE1 - SEQUENTIAL FILE CONTAINING SITE II INPUT 
DECK, EACH DESCRIBING 150 POINTS OR LESS, 
WITH INPUT DECKS SEPARATED BY E-O-F MARKS. 
SEE SUBROUTINE BLDARR. 

SCRATCH: 

TAPE4 - SEQUENTIAL FILE CONTAINING CARD IMAGES OF 
SITE II INPUT DECKS BEFORE PROCESSING BY 
BLDARR. 
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SUBROUTINE CHKREC 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OBJECTIVE: 

TO DETERMINE IF A MESSAGE, REGARDING USER-DEFINED PRO- 
FILES, IS WRITTEN ON THE INPUT FILE (TAPE5) 


METHOD: 

SUBROUTINE CHKREC EVALUATES TAPE5 TO DETERMINE THE 
NUMBER OF RECORDS BETWEEN THE CARD IMAGE CONTAINING 
THE NOISE METRIC VALUE AND THE BEGINNING OF THE 
XY-DATA POINTS 


PARAMETERS: 

INPUT: 

NONE 

OUTPUT: 

ICASE: 8-WORD ARRAY CONTAINING CASE NAME IN WORDS 5 
THROUGH 8 

ISW: INTEGER VARIABLE USED TO INDICATE THE NUMBER 

OF RECORDS TO BE SKIPPED PRIOR TO THE READING 
OF THE XY DATA POINTS 


CALLING ROUTINE: 

AUPHS1 (MAIN PROGRAM) 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

IBLANK: INTEGER ARRAY USED TO CONTAIN THE POTENTIAL 
SPECIAL MESSAGE 

J: INTEGER VARIABLE USED AS LOOP COUNTER 

KEOF: INTEGER VARIABLE USED AS A FLAG TO INDICATE AN 

EARLY END-OF- INPUT 
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FILES USED: 

INPUT: 

TAPE5 - SEQUENTIAL FILE CONSISTING OF VARIABLE- 

LENGTH RECORDS; SUBROUTINE CHKREC IS CONCERNED 
WITH THE FORMAT OF THE FIFTH Z-TYPE RECORD ON 
TAPE5; WITH NO SPECIAL MESSAGE, THIS RECORD 
IS BLANK; WITH SPECIAL MESSAGE, RECORD CON- 
SISTS OF A SPECIAL MESSAGE, FOLLOWED BY A 
VARIABLE NUMBER OF EMPTY RECORDS 

OUTPUT: 

NONE 

SCRATCH: 

NONE 
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SUBROUTINE RDCONTR 


OBJECTIVE: 

THIS ROUTINE READS THE HEADER DATA AND THE X, Y POINTS FOR 
A NOISE CONTOUR AS OUTPUT BY INM. 


PARAMETERS: 

INPUT: 

ISW: INTEGER VARIABLE USED TO POSITION CORRECTLY 

TO XY DATA POINTS ON TAPE5 

OUTPUT: 

LINE: STORAGE FOR NOISE METRIC AND VALUES. 

X: X POINTS FOR CONTOUR 

Y: Y POINTS FOR CONTOUR 

NPTS: THE NUMBER OF X,Y POINTS. 

KEOF: SET TO 1 TO INDICATE EOF ON INPUT FILE. 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

EOF : SYSTEM END-OF'-FILE CHECK. 


LOCAL VARIABLE DICTIONARY: 

A, B, C, D 
E, F,I,J 

K,L,M : DUMMY VARIABLES FOR I/O LIST ONLY. 

IJ : IMPLIED DO-LOOP INDEX. 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

TAPE5 - SEQUENTIAL INPUT FILE CONTAINING THE CONTOUR X f Y 
POINTS AND HEADER DATA. (SEE MAIN PROGRAM) 

OUTPUT: 

NONE 
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SUBROUTINE INIT 
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OBJECTIVE: 

THIS ROUTINE READS THE SUPPLEMENTAL DATA FILE, TAPE2. 


PARAMETERS: 

INPUT: 

NONE 

OUTPUT: 

AIRPORT CENTROID LATITUDE. 
AIRPORT 1 CENTROID LONGITUDE. 
NORTH/SOUTH DISPLACEMENT. 
EAST/WEST DISPLACEMENT. 

X AXIS TRANSLATION DISTANCE. 
Y AXIS TRANSLATION DISTANCE. 
PLOT AXIS LENGTH. 

AIRPORT TITLE. 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 
NONE 


LONG 

NSDIS 

EWDIS 

H 

K 

FACT 

ITITLE 


FUNCTIONS USED: 

EOF : SYSTEM END-OF-FILE TEST. 


LOCAL VARIABLE DICTIONARY: 
NO LOCAL VARIABLES 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

TAPE2 - SEQUENTIAL INPUT FILE IN DISPLAY CODE FORMAT 
SEE MAIN PROGRAM FOR CARD BY CARD FORMAT. 
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SUBROUTINE SITE2 
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OBJECTIVE: 

THIS ROUTINE DRIVES THE POINT REDUCTION SOFTWARE. ITS 
OBJECTIVE IS TO DELETE AS MANY POINTS AS POSSIBLE FROM 
A NOISE CONTOUR WITHOUT CHANGING THE TOTAL AREA THAT THE 
CONTOUR ENCLOSES BY MORE THAN 1 PERCENT. THE MINIMUM 
OBJECTIVE IS TO REDUCE THE CONTOUR TO LESS THAN 150 PI'S. 

METHOD: 

THE ROUTINE SETS THE REDUCTION TARGET TO 65 POINTS AND 
CALCULATES THE AREA OF THE ORIGINAL CONTOUR. THE POINT 
REDUCTION ROUTINE IS CALLED ITERATIVELY UNTIL THE AREA 
RESTRICTION IS MET. THE TARGET NUMBER OF POINTS IS 
INCREASED BY 5 ON EACH ITERATION. 


PARAMETERS : 

INPUT: 

X : THE ORIGINAL CONTOUR X POINTS. 

Y : THE ORIGINAL CONTOUR Y POINTS. 

N : THE NUMBER OF POINTS IN THE ORIG. CONTOUR. 

OUTPUT: 

XNEW : THE X POINTS AFTER REDUCTION. 

YNEW : THE Y POINTS AFTER REDUCTION. 

NEWNUM : THE NUMBER OF POINTS AFTER POINT REDUCT. 


CALLING ROUTINE: 
PHASE1 


SUBPROGRAMS USED: 

NEWCRV : POINT REDUCTION ALGORITHM. 

FUNCTIONS USED: 

AREA : CALCULATES THE AREA WITHIN A CLOSED POLYGON. 

ABS : ABSOLUTE VALUE FOR REAL NUMBERS. 

LOCAL VARIABLE DICTIONARY: 

AREA1 : THE ORIGINAL CONTOUR AREA. 

AREA2 : THE AREA FOR THE CONTOUR AFTER POINT REDUCTION. 

DIF : THE DIFFERENCE IN AREA FOR THE ORIGINAL AND 

POINT REDUCED CONTOURS. 

INT : THE NUMBER OF POINTS TO RELAX POINT REDUCTION 

BY TO MEET THE 1 PERCENT CHANGE IN AREA LIMIT. 
ITER : THE NUMBER OF ITERATIONS THAT SUBROUTINE NEWCRV 

PERFORMED TO DELETE POINTS. 

PRCNT : THE CALCULATED PERCENT CHANGE IN AREA AFTER 

POINT REDUCTION. 
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SUBROUTINE NEWCRV 
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OBJECTIVE: 

TO DELETE THE MOST UNNECESSARY POINTS BELOW A GIVEN 
MAXIMUM ON A NOISE CONTOUR. POINTS DELETED WILL BE THOSE 
WHICH TEND TO BE ON STRAIGHT LINE SEGMENTS OR BROAD 
CURVES. 


METHOD: 

OE 3 CONTIGUOUS POINTS A, B AND C, POINT B WILL BE DELETED 
WHEN THE ABSOLUTE VALUE OF THE DIFFERENCE IN SLOPES FOR 
LINES AB AND BC IS LESS THAN A GIVEN LIMIT. THE SIGNED 
DIFFERENCES ARE ACCUMULATED SUCH THAT WHEN THE MAGNITUDE 
OF THE ACCUMULATION BECOMES GREATER THAN THE LIMIT THE 
POINT TO BE DELETED IS SAVED. THIS PROHIBITS EXCESSIVE 
CONTOUR DISTORTION OVER BROAD CURVES IN ONE DIRECTION. 

IF ENOUGH POINTS CANNOT BE DELETED AFTER PROCESSING THE 
CONTOUR THE LIMIT IS RELAXED AND A NEW ITERATION IS 
BEGUN. AS THE LIMIT IS RELAXED MORE POINTS WILL BE 
DELETED AND MORE DISTORTION OF THE CONTOUR WILL BE 
INTRODUCED. 


PARAMETERS: 

INPUT: 

X : 

Y : 

N : 

MAXNUM : 

OUTPUT: 

XNEw : 
YNEW : 
NEWNUM : 
ITER : 


X POINTS OF THE CONTOUR TO BE REDUCED. 

Y POINTS OF THE CONTOUR TO BE REDUCED. 

THE NUMBER OF X,Y POINTS IN CONTOUR TO BE 
REDUCED. 

THE MAXIMUM NUMBER OF POINTS TO BE LEFT IN 
THE REDUCED CONTOUR. 

THE REDUCED X POINTS. 

THE REDUCED Y POINTS. 

THE NUMBER OF X,Y POINTS AFTER REDUCTION. 
THE NUMBER OF SEPARATE COMPLETE LOOPS 
AROUND THE CONTOUR. IE, THE NUMBER OF 
RELAXATIONS OF THE SLOPE TOLERANCE. 


CALLING ROUTINE: 
SITE2 


FUNCTIONS USED: 

IOCTNT : USED TO DETERMINE WHAT COMPASS OCTANT AN X,Y 

POINT IS IN. 
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LOCAL VARIABLE DICTIONARY: . 

B : ORDINAL POINTER TO THE 1ST OF THE 3 POINTS USED 

TO CALCULATE THE SLOPES. 

DELPTS : ARRAY STORAGE FOR ORDINALS OF DELETED POINTS. 

DIF : THE DIFFERENCE IN THE SLOPES, BT AND TTP1. 

I : DO LOOP INDEX OVER THE TOTAL PT. ORDINALS FOR 

THE ORIGINAL CONTOUR. 

INC : THE RELAXATION INCREMENT VALUE FOR THE SLOPE 

DIFFERENCE LIMIT. 

KI : TOE OCTANT NUMBER OF THE MIDPOINT OF THE 3 

POINTS BEING EVALUATED. (B,T,TPl) 

KJ : TOE OCTANT NUMBER OF THE 3RD POINT, TP1 OF THE 

3 POINTS BEING EVALUATED. 

KPT : ORDINAL POINTER FOR THE ARRAY OF NEW POINTS 

AFTER REDUCTION. 

LAST : TOE LAST ORDINAL ON THE ORIGINAL CONTOUR. 

MAXCHG : THE LIMIT FOR DIFFERENCE IN SLOPES WHEN A POINT 

IS OR IS NOT DELETED. 

Ml : THE SLOPE OF LINE SEGMENT BT. 

M2 : THE SLOPE OF LINE SEGMENT TTP1.(TP1= T + 1) 

NEXT : NEXT AVAILABLE SLOT PTR FOR DELETED POINTS 

ARRAY, DELPTS. 

NPT : NEXT AVAILABLE SLOT PTR FOR NEW POINTS ARRAY. 

SUM : ACCUMULATOR FOR SIGNED SLOPE DIFFERENCES . 

T : THE MIDDLE POINT ORDINAL FOR THE 3 POINTS UNDER 

SCRUTINY. T IS THE DELETION POINT. 

TP1 : THE 3RD POINT ORDINAL, T + 1. 

XZ : DELTA X FOR THE BT SLOPE. 

XZ1 : DELTA X FOR THE TIPI SLOPE. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE NTHSTH 
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OBJECTIVE: 

TO CONVERT X,Y POINTS FROM FEET TO STATUTE MILES IN A 
SITE II COMPATIBLE FORMAT USING NORTH/SOUTH/EAST/WEST 
DESIGNATORS AND TO OUTPUT THE CONVERTED POINTS AS SITE II 
"D" CARDS. 


METHOD: 

THE POINTS ARE CONVERTED IROM FEET TO STATUTE MILES BY 
DIVIDING THE FEET UNITS BY 5280. Y POINTS THAT ARE 
POSITIVE ARE GIVEN AN "N" (NORTH) INDICATOR; NEGATIVE Y'S 
ARE GIVEN AN "S" (SOUTH) INDICATOR. THE SAME LOGIC 
APPLIES TO X POINTS USING E/W (EAST/WEST) INDICATORS. 


PARAMETERS: 

INPUT: 

X : THE X PTS TO BE CONVERTED. 

Y : THE Y PTS TO BE CONVERTED. 

N : THE NUMBER OF X,Y POINTS. 


CALLING ROUTINE: 
POLY 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

ABS : ABSOLUTE VALUE FOR REAL NUMBERS. 
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LOCAL VARIABLE DICTIONARY: 

I : DO LOOP INDEX TO COLUMNS OF ARRAY JBUF, A 3 

CO-ORDINATE BUFFER (1"D" CARD COMPLEMENT) . 

J : IMPLIED DO LOOP INDEX FOR COLUMNS OF JBUF. 

JBUF : BUFFER ARRAY FOR 3 CO-ORDINATES FOR A "D" CARD. 

K : IMPLIED DO LOOP INDEX FOR ROWS OF JBUF WHEN 

THERE ARE LESS THAN 3 POINTS (LAST "D" CARD) . 
KPT : TOTAL CO-ORDINATE COUNTER TO PREVENT WRITING 

: ONE OR TOO EXTRA POINTS ON LAST "D" CARDS. 

KTRI : FORMAT REPETITION COUNT. 

L : IMPLIED DO LOOP INDEX FOR COLUMNS OF JBUF 

XMILES : THE ABSOLUTE VALUE IN MILES OF THE CURRENT 

X POINT. 

YMILES : SAME AS XMILES BUT FOR Y POINTS. 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

NONE 

SCRATCH: 

TAPE4 - SEQUENTIAL FILE FOR SITE II INPUT DECKS. 

IT IS SUBSEQUENTLY REPROCESSED FOR 150 POINT 
RESTRICTION VIOLATIONS. 
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SUBROUTINE TKNSLT 


OBJECTIVE: 

TO TRANSLATE TOE AXIS FOR TOE NOISE CONTOURS. 


METHOD: 

THE TRANSLATION VALUES H AND K ARE SUBTRACTED FROM EACH 
X AND Y POINT RESPECTIVELY. 


PARAMETERS : 

INPUT: 


X 

Y 

NPTS 

H 

K 


X PI’S FOR CONTOUR. 

Y ITS FOR CONTOUR. 

TOE NUMBER OF X,Y PTS. 

THE X TRANSLATION E ACTOR. 
THE Y TRANSLATION FACTOR. 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

I : DO LOOP INDEX OVER THE NUMBER OF PTS. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 


A- 15 


T 



FUNCTION AREA 
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OBJECTIVE: 

TO CALCULATE THE AREA IN SQ UNITS OF A CLOSED POLYGON. 


METHOD: 

A SPECIAL BSRM OF GREEN'S THEOREM IS EMPLOYED: 
N 

AREA = .5 * SUM (X *Y - X *Y ) 

1=0 I 1+1 1+1 I 

WHERE X = X Y = Y 

N+l 0 N+l 0 


PARAMETERS: 

INPUT: 

X : THE CONTOUR X POINTS. 

Y : THE CONTOUR Y POINTS. 

N : THE NUMBER OF X,Y POINTS. 

OUTPUT: 

AREA : FUNCTION VALUE RETURNED FOR THE CONTOUR 

AREA. 


CALLING ROUTINE: 
SITE2 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

ABS : ABSOLUTE VALUE. 


LOCAL VARIABLE DICTIONARY: 

AREA : SQUARE UNITS INSIDE CONTOUR. 

I : DO LOOP INDEX OVER RANGE OF X,Y POINTS. 

K : THE 1+1 ORDINAL. 

SUM : REAL VARIABLE USED IN COMPUTING CONTOUR 

AREA 
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SUBROUTINE OCTANT 
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OBJECTIVE: 

THIS ROUTINE CO-ORDINATES THE POLYGON DEFINITIONS. 


PARAMETERS: 

INPUT: 

X : THE CONTOUR X POINTS. 

Y : THE CONTOUR Y POINTS. 

NEWPTS : THE NUMBER OF X, Y POINTS AFTER POINT DEL. 

LDN : THE NOISE METRIC OF THE CONTOUR. (LDN) 

MET : TOE NOISE METRIC VALUE. (55 - 85) 

INDEX : THE CONTOUR COUNTER. 

SCALE : THE PLOTTING SCALE FACTOR. 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 

OCTTAB : BUILDS A TAbLE OF START AND END POINTS OF LINE 

SEGMENTS FOR EACH OCTANT. 

POLY : BUILDS THE POLYGONS USING THE TABLE BUILT BY 

OCTTAB. 

ENRICH : CHECKS FOR OCTANTS SPANNED WITH NO X,Y POINTS 

AND PUTS IN POINTS WHEN NEEDED. 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

ITABLE : THE START POINT/END POINT TABLE FOR LINE 

SEGMENTS IN EACH OCTANT. 

KCHANGE: FLAG THAT INDICATES POINTS WERE PUT IN AND 

ITABLE MUST BE REDEFINED BY OCTTAB ROUTINE. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE ENRICH 



OBJECTIVE: 

SUBROUTINE ENRICH IS USED TO EVALUATE CERTAIN POINTS TO 
DETERMINE IF ADJACENT POINTS ARE IN EITHER (1) THE SAME 
OCTANT OR (2) ADJACENT OCTANTS; IF 2 ADJACENT POINTS NOT 
IN THE SAME OCTANT ARE NOT IN ADJACENT OCTANTS, ROUTINE 
ADDPTS IS INVOKED TO ADD POINTS IN ORDER TO DEAL WITH 
THE SITUATION 


PARAMETERS: 

INPUT: 

ITABLE : INTEGER ARRAY USED TO STORE THE START/END 
POINTS FOR LINE SEGMENTS IN EACH OCTANT 
NEWPTS : THE NUMBER OF XY POINTS AFTER POINT 
DELETION 

X : REAL ARRAY CONTAINING THE X-COORDINATES 

OF CONTOUR POINTS PRIOR TO POSSIBLE POINT 
ENRICHMENT 

Y : REAL ARRAY CONTAINING THE Y-COORDINATES 

OF CONTOUR POINTS PRIOR TO POSSIBLE 
ENRICHMENT 

OUTPUT: 

KCHANGE : INTEGER VARIABLE USED AS FLAG TO INDICATE 
THAT POINTS WERE PUT IN AND ITABLE MUST 
BE REDEFINED BY OCTTAB ROUTINE 
NEWPTS : INTEGER VARIABLE REPRESENTING THE NUMBER 
OF XY POINTS AFTER POINT ENRICHMENT 
X : REAL ARRAY CONTAINING X-COORDINATES OF 

CONTOUR POINTS AFTER POSSIBLE ENRICHMENT 

Y : REAL ARRAY CONTAINING Y-COORDINATES OF 

CONTOUR POINTS AFTER POSSIBLE ENRICHMENT 


CALLING ROUTINE: 

OCTANT 

SUBPROGRAMS USED: 

ADDPTS: ROUTINE USED TO ADD POINTS ON A CONTOUR WHEN A 
SPECIFIED CONDITION, AS DESCRIBED IN "OBJECTIVE" 
SECTION 

SYSTEM: STANDARD FTN SUBROUTINE USED TO TERMINATE 

PROGRAM EXECUTION AND WRITE A DAYFILE MESSAGE 
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FUNCTIONS USED: 

INTABLE : FUNCTION USED TO DETERMINE THE OCTANT OF AN 
ORDINAL, IF THE ORDINAL IS IN THE START/END 
TABLE 

IOCTNT : FUNCTION USED TO DETERMINE WHAT OCTANT AN 
X, Y POINT RESIDES IN 

LOCAL VARIABLE DICTIONARY: 

A : REAL ARRAY USED FOR TEMPORARY STORAGE OF THE 

X-COORDINATES OF CONTOUR POINTS 

B : REAL ARRAY USED FOR TEMPORARY STORAGE OF THE 

Y-COORDINATES OF CONTOUR POINTS 

I : INTEGER VARIABLE USED AS LOOP COUNTER AND ARRAY 

INDEX 

IOCT : INTEGER VARIABLE USED TO INDICATE IF ORDINAL 
IS IN THE START/END TABLE 

IOCTM1 : INTEGER VARIABLE REPRESENTING AN ADJACENT 
OCTANT (COUNTERCLOCKWISE) 

IOCTP1 : INTEGER VARIABLE REPRESENTING AN ADJACENT 
OCTANT (CLOCKWISE) 

IPTR : INTEGER USED TO COUNT POINTS IN CONTOUR 

J : INTEGER VARIABLE USED TO REPRESENT THE ORDINAL 

OF THE NEXT CLOCKWISE POINT 

KOCT : INTEGER VARIABLE REPRESENTING THE OCTANT OF 
THE ADJACENT POINT 

NPTS : 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

NONE 

SCRATCH: 

NONE 


A- 19 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION INTABLE 


OBJECTIVE: 

FUNCTION INTABLE IS USED TO DETERMINE THE OCTANT OF AN 
ORDINAL 


PARAMETERS: 

INPUT: 

KORD : ORDINAL OF XY POINT 

ITABLE : THE START/END POINT TABLE FOR LINE SEGMENTS 
IN EACH OCTANT 

OUTPUT: 

NONE 


CALLING ROUTINE: 
ENRICH 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

I : INTEGER VARIABLE .USED AS LOOP COUNTER AND ARRAY INDEX 
J : INTEGER VARIABLE USED AS LOOP COUNTER AND ARRAY INDEX 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

NONE 

SCRATCH: 

NONE 
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SUBROUTINE ADOPTS 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OBJECTIVE : 

SUBROUTINE USED TO ADD POINTS ON A CONTOUR WHEN A 
SPECIFIC CONDITION, AS DESCRIBED IN ENRICH, IS ENCOUNTERED 


PARAMETERS: 

INPUT: 

A : REAL ARRAY CONTAINING X VALUES OF CONTOUR 
POINTS 

B : REAL ARRAY CONTAINING Y VALUES OF CONTOUR 
POINTS 

IPTR : INTEGER USED TO COUNT POINTS IN A CONTOUR 

NPTS : INTEGER SET TO VALUE OF 2 

XST : X-COORDINATE OF POINT PRECEDING ANOTHER 

POINT THAT IS NOT IN THE SAME, OR ADJACENT, 
OCTANT 

YST : Y-COORDINATE OF POINT PRECEDING ANOTHER 

POINT THAT IS NOT IN THE SAME, OR ADJACENT, 
OCTANTS 

XEND : X-COORDINATE OF POINT THAT HAS A PRECEDING 
POINT THAT SPANS MORE THN 1 OCTANT 

YEND : Y-COORDINATE OF POINT THAT HAS A PRECEDING 
POINT THAT SPANS MORE THAN 1 OCTANT 

OUTPUT: 

A : REAL ARRAY CONTAINING X VALUES OF CONTOUR 
POINTS AFTER POINT ENRICHMENT 

B : REAL ARRAY CONTAINING Y VALUES OF CONTOUR 
POINTS AFTER POINT ENRICHMENT 


CALLING ROUTINE: 
ENRICH 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

ABS : STANDARD FTN FUNCTION USED TO CALCULATE ABSOLUTE 
VALUE 

FLOAT : STANDARD FTN FUNCTION USED TO CONVERT INTEGER 
TO REAL NUMBER 
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c 

c 
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LOCAL VARIABLE DICTIONARY: 

DX : DIFFERENCE IN THE VALUES OF THE X-COORDINATES 
OF THE POINTS UNDER EVALUATION 
DY : DIFFERENCE IN THE VALUES OF THE Y-COORDINATES 
OF THE POINTS UNDER EVALUATION 
I : INTEGER VARIABLE USED AS LOOP COUNTER 
X : X-COORDINATE OF NEW CONTOUR POINT FOUND 
XINT : REAL VARIABLE REPRESENTING INCREMENT ON THE 
STARTING X VALUE, WHEN COMPUTING NEW VALUE 
Y : Y-COORDINATE OF NEW CONTOUR POINT FOUND 
YINT : REAL VARIABLE REPRESENTING INCREMENT OF 

STARTING X-VALUE WHEN COMPUTING NEW POINT 
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SUBROUTINE OCTTAB 
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c 
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c 
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c 
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c 

c 

c 
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OBJECTIVE: 

THIS ROUTINE BUILDS A TABLE OF START AND END POINTS FOR 
EACH LINE SEGMENT IN EACH OCTANT FOR THE SUBSEQUENT 
POLYGON DEFINITIONS . 


METHOD: 

THE CONTOUR POINTS ARE ARRANGED IN A CLOCKWISE FASHION IN 
ASCENDING ORDINAL ORDER. IT IS TYPICAL THAT THE 1ST 
POINT FALLS IN THE MIDDLE OF AN OCTANT. THIS PROCESSOR 
LOOKS BACKWARD (COUNTERCLOCKWISE) FROM THE 1ST POINT 
UNTIL IT FINDS A POINT OUTSIDE THE OCTANT THAT THE 1ST 
POINT OCCUPIES. THE LAST COUNTERCLOCKWISE POINT THAT IS 
INSIDE THE STARTING OCTANT IS INSERTED IN THE TABLE AS 
INDEXED BY THE OCTANT (1-8) . THE COLUMNS OF THE TABLE 
ARE ARRANGED AS ODD FOR START POINTS AND EVEN FOR END 
POINTS. FOR EXAMPLE, THE SECOND LINE SEGMENT IN OCTANT 4 
WOULD BE STORED IN ROW 4, WITH THE START POINT IN COLUMN 3 
AND THE END POINT IN COLUMN 4. THE ARRAY NXTCOL IS USED TO 
POINT TO THE NEXT AVAILABLE COLUMN ORDINAL IN EACH OCTANT 
OF THE TABLE. HAVING FOUND THE REAL STARTING POINT BY 
BACKING UP AS DESCRIBED ABOVE, THE END POINT IS FOUND IN 
THE SAME MANNER BY MOVING AHEAD (CLOCKWISE) UNTIL THE END 
POINT IS FOUND AND STORED IN THE TABLE. AT THIS TIME THE 
CURRENT OCTANT BECOMES THE OCTANT THAT THE 1ST OUTSIDE 
POINT WAS FOUND WITHIN. THE POINT IS STORED IN THE TABLE 
AS A START POINT AND THE CLOCKWISE PROCESS IS BEGUN AGAIN 
BY MOVING CLOCKWISE UNTIL THE NEXT OUTSIDE POINT IS 
FOUND. THIS PROCESSING IS CONTINUED UNTIL THE ORIGINAL 
STARTING POINT IS ENCOUNTERED (NOT NECESSARILY ORDINAL 1) 
INDICATING 1 COMPLETE CLOCKWISE TRIP AROUND THE CONTOUR. 


PARAMETERS: 

INPUT: 

X : THE X POINTS FOR THE POINT REDUCED CONTOUR 

Y : THE Y POINTS FOR THE POINT REDUCED CONTOUR 

NPTS : THE NUMBER OF X,Y POINTS. NOTE THAT POINT 

1 = POINT NPTS. 

OUTPUT: 

TABLE : THE STARTPOINT/ENDPOINT LINE SEGMENT, 

OCTANT TABLE. 
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c 

c 

c 

c 

c 
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CALLING ROUTINE: 
OCTANT 


SUBPROGRAMS USED: 

NONE 

FUNCTIONS USED: 

IOCTNT : DETERMINES THE OCTANT WHERE AN X,Y POINT 

RESIDES. 


LOCAL VARIABLE DICTIONARY: 

BEGOCT : THE ORDINAL OF THE BEGINNING OCTANT (WHERE FT 

NUMBER 1 RESIDES) . 

CUROCT : THE ORDINAL OF THE OCTANT WHERE PTS ARE BEING 

PROCESSED TO FIND AN OUTSIDE PT. 

I : DO LOOP INDEX 

INDEX : THE CONTOUR ORDINAL POINTER FOR CLOCKWISE AND 

COUNTERCLOCKWISE PROCESSING. 

LSTPT : THE ORDINAL OF THE LAST POINT IN A LINE SEGMENT 

WHICH IS IN THE SAME OCTANT AS THE OTHER POINTS 

NXTCOL : THE VECTOR OF PTRS TO THE NEXT AVAILABLE PAIR 

OF COLUMNS FOR A ST/END POINT PAIR FOR THE ITH 
OCTANT. 

NXTOCT : THE OCTANT ORDINAL FOR THE NEXT CLOCKWISE PT ON 

A LINE SEGMENT WHICH IS COMPARED WITH THE 
VARIABLE CUROCT. 

STRTPT : THE TRUE STARTING POINT OF THE 1ST LINE SEGMENT 

FOUND BY BACKING UP FROM PT NUMBER 1. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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FUNCTION IOCTNT 
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OBJECTIVE: 

TO DETERMINE WHAT OCTANT AN X,Y POINT RESIDES IN. 


METHOD: 

THE OCTANTS ARE ORDERED CLOCKWISE AS FOLLOWS USING THE 
STANDARD CARTESIAN QUANDRANTS: 

OCTANT 1 THE UPPER HALF OF QUADRANT IV 

2 THE LOWER HALF OF QUADRANT IV 

3 THE LOWER HALF OF QUADRANT III 

4 THE UPPER HALF OF QUADRANT III 

5 THE DOWER HALF OF QUADRANT II 

6 THE UPPER HALF OF QUADRANT II 

7 THE UPPER HALF OF QUADRANT I 

8 THE LOWER HALF OF QUADRANT I 

SPECIAL CASES ARE TESTED FOR (X = Y OR X = Y = 0) FIRST 
THEN THE QUADRANT IS TESTED FOR. ONCE THE QUADRANT IS 
DETERMINED, IT IS DETERMINED ON WHAT SIDE OF THE LINE 
Y = X OR Y = -X THE POINT FALLS. THE SPECIAL CASES WHERE 
THE POINT FALLS ON AN OCTANT LINE OR ON THE ORIGIN CAUSE 
THE OCTANT TO BE DEFINED AS THE LAST OCTANT VALUE. 


PARAMETERS: 

INPUT: 

X : THE X POINTS OF THE CONTOUR. 

Y : THE Y POINTS OF THE CONTOUR. 

N : THE NUMBER OF POINTS IN THE CONTOUR. 

INDEX : THE ORDINAL OF THE POINT TO BE EXAMINED. 

LSTOCT : THE CURRENT OCTANT ORDINAL. 

OUTPUT: 

NONE 


CALLING ROUTINE: 
NEWCRV 
OCTTAB 
POLY 
DSTNCE 
INTSCT 


FUNCTIONS USED: 

ABS : ABSOLUTE VALUE ROUTINE. 
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SUBROUTINE POLY 
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OBJECTIVE: 

THIS ROUTINE FORMS POLYGONS FROM THE START/END POINT 
TABLE DEFINED BY SUBROUTINE OCTTAB. THE POLYGONS ARE 
FORMED BY THE CONTOUR LINES AND THEIR INTERSECTION WITH 
THE OCTANT LINES. WHEN A POLYGON IS DEFINED IT IS OUTPUT 
AS A SITE II "C" CARD FOLLOWED BY THE REQUIRED NUMBER OF 
"D" CARDS. LASTLY THE POLYGON IS PLOTTED IF THE CONTOUR 
COUNTER IS ODD. 


METHOD: 

THE TABLE OF START AND END POINTS IS 1ST PROCESSED TO 
COMPUTE THE DISTANCE FROM THE ORIGIN TO EACH START AND 
END POINT ON AN OCTANT LINE. TO DO THIS THE EQUATION OF 
A LINE IS DEFINED USING A START/END POINT AND ITS 
NEIGHBORING OUTSIDE POINT. THIS LINE MUST INTERSECT AN 
OCTANT LINE AND IT IS THIS INTERSECTION POINT THAT IS 
USED FOR DISTANCE TO ORIGIN CALCULATIONS. THE DISTANCES 
ARE STORED IN THE TABLE DSTS. ANOTHER TABLE, IUSED, IS 
USED TO MARK A START/END POINT AS USED OR UNUSED SINCE 
THE ALGORITHM WILL DO THINGS LIKE SEARCH AWAY FROM THE 
ORIGIN FOR THE 1ST UNUSED START PT OR TOWARDS THE ORIGIN 
FOR THE 1ST UNUSED END PT. THE DISTANCES WERE CALCULATED 
TO AFFORD THIS CAPABILITY. THE POLYGON GENERATION 
ALGORITHM CAN BE STATED AS THE FOLLOWING 4 RULES: 

NOTE: THE CURRENT OCTANT LINE IS THE OCTANT LINE WHOSE 
ORDINAL = THE OCTANT NUMBER BEING PROCESSED. THE 
POSITIVE X AXIS IS ORDINAL 1. THE OTHER OCTANT 
LINES ARE NUMBERED ASCENDINGLY FROM NO.l IN A 
CLOCKWISE DIRECTION. 

1. START PROCESSING AN OCTANT BY TAKING POINTS STARTING 
AT THE 1ST START POINT IN THE TABLE UNTIL AN END 
POINT IS ENCOUNTERED. 

2. WHEN AN ENDPOINT IS ENCOUNTERED ON THE CURRENT OCTANT 
LINE FIND THE NEXT START POINT MOVING AWAY FROM THE 
ORIGIN (A START PT WHOSE DISTANCE IS GREATER THAN THE 
END PT JUST ENCOUNTERED AND WHOSE DISTANCE IS THE 
MINIMUM OF OTHER QUALIFYING START FTS) . IF THE 
QUALIFYING START FT HAS ALREADY BEEN USED WE HAVE 
CLOSURE (A COMPLETE POLYGON) . IF IT HAS NOT BEEN USED 
IT IS SELECTED AS A MEMBER OF THE CURRENT POLYGON AND 
THE "USED" TABLE IS UPDATED. 
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3. WHEN AN ENDPOINT IS ENCOUNTERED ON THE NEXT OCTANT 
LINE FIND THE NEXT START POINT MOVING TOWARD THE 
ORIGIN (A START PT WHOSE DISTANCE IS LESS THAN THE END 
FT JUST ENCOUNTERED AND WHOSE DISTANCE IS THE MAXIMUM 
OF OTHER QUALIFYING START PTS) . IF A START POINT CAN 
NOT BE FOUND THEN THE ORIGIN IS TAKEN AS THE NEXT PT. 
AND RULE NUMBER 2 IS APPLIED. 

4. WHEN A START/END POINT IS SELECTED IT IS NO LONGER 
AVAILABLE FOR LATER SELECTION. 


c 

PARAMETERS: 



c 

INPUT: 



c 

TABLE 

THE 

START/END POINT OCTANT TABLE. 

c 

X 

THE 

X POINTS OF THE CONTOUR. 

c 

Y 

THE 

Y POINTS OF THE CONTOUR. 

c 

N 

THE 

NUMBER OF X,Y POINTS. 

c 

MET 

THE 

TABLE OF LDN METRICS. 

c 

LDNPT 

pro 

TO CURRENT LDN METRIC. 

c 

LDN 

THE 

METRIC NAME. ("LDN") . 

c 

c 

SCALE 

THE 

PLOTTING SCALE FACTOR. 

c 

c 

i 

CALLING ROUTINE: 



c 

c 

OCTANT 



c 

c 

SUBPROGRAMS USED: 



c 

D6TNCE : CALCULATES THE DISTANCE FROM THE ORIGIN 


INTSCT 


NTHSTH 


PLTPOL 


FUNCTIONS USED: 
AREA 
IOCTNT 


START/END POINTS. 

FINDS THE INTERSECTION X, Y POINT OF AN OCTANT 
LINE AND A LINE DEFINED BY 2 X,Y POINTS WHICH 
STRADDLE THE OCTANT LINE. 

CONVERTS THE POLYGON POINTS TO SITE II 
NORTH/ SOUTH/EAST/WEST DISPLACEMENTS IN MILES. 
ROUTINE WHICH PLOTS POLYGON X,Y POINTS. 


CALCULATES THE AREA OF A POLYGON. 

DETERMINES THE OCTANT NUMBER FOR A GIVEN X,Y 
POINT. 
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c 

c 

c 

c 
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LOCAL VARIABLE DICTIONARY: 

A : THE POLYGON AREA IN SQUARE MILES, OUTPUT AS F6.1 

ON THE SITE II "C" CARD. 

DELTAX : WHEN AN INTERSECTION PT IS CALCULATED ON AN 

OCTANT LINE, ITS PROXIMITY TO THE ALREADY 
EXISTING NEIGHBOR PT IS CALCULATED IN TERMS OF 
DELTAX AND DELTAY (THE DIFFERENCE IN X'S AND 
Y'S) . 

DELTAY : SEE DELTAX. 

DSTS : A TABLE OF START/END POINT DISTANCES FROM THE 

ORIGIN. 

ENDDIS : A DISTANCE TO A SELECTED END POINT (FROM THE 


ORIGIN) . 

I : DO LOOP VARIABLE FOR INITIALIZATION. 

ICOL : COLUMN INDEX TO LINE SEGMENT ST/END POINT TABLE 

ICOUNT : COUNTER FOR TOE NUMBER OF TIMES THIS ROUTINE IS 

CALLED. USED TO PLOT POLYGONS ON ODD COUNTS. 

IF LAG : SET TO 1 TO INDICATE THAT THIS IS AN ODD COUNT 

CALL. 

INPT : THE "INSIDE POINT" WHEN 2 POINTS STRADDLE AN 

OCTANT LINE. THE INSIDE POINT IS THE POINT 
WITHIN THE OCTANT BEING ANALYSED. 

IOCT : THE ARRAY INDEX CORRESPONDING TO OCTANTS. 

IPTR : THE NEXT AVAILABLE SLOT FOR X, Y POINT STORAGE 

AS THE POLYGON IS BEING FORMED. 

IRES : THE RESULT OF DIVIDING ICOUNT BY 2. WHEN 

ICOUNT IS ODD THE REMAINDER IS TRUNCATED 
SUCH THAT 2* IRES WILL NOT = ICOUNT. 

I USED : AN ARRAY USED TO MARK ST/END POINTS AS USED OR 

UNUSED DURING POLYGON DEFINITION. 

J : MISCELLANEOUS DO LOOP COUNTER. 

JPTR : THE NEXT AVAILABLE SLOT FOR A VECTOR OF 

QUALIFYING ST/END POINTS TO BE COMPARED FOR 
DISTANCES FROM THE ORIGIN. 

KOCT : THE OCTANT ORDINAL OF AN INTERSECTION POINT ON 

AN OCTANT LINE. 

KPT : WORKING POINTER TO THE ST/END POINTS COLUMNS. 

KVEK : VECTOR FOR STORAGE OF ST/END POINTS TO BE 

COMPARED FOR MINS AND MAXES. 

LAST : LAST POINT (END POINT) IN A LINE SEGMENT. 

MVEK : SAME AS KVEK. 

NPTR : WORKING POINTER TO EITHER AN INTERSECTION PT OR 

AN EXISTING POINT EXTREMELY CLOSE TO THE 
INTERSECTION POINT. 

NUM : ON AN EVEN VALUE OF ICOUNT NUM WILL BE EQUAL TO 

ICOUNT. ON ODD VALUE NUM WILL NOT EQUAL ICOUNT 
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c 

OUTPT 

: 

SEE INPT. THE OUTSIDE POINT. 


c 

WHERE 

; 

AN INTEGER ARRAY OF OCTANT NUMBERS THAT 


c 



CORRESPOND TO THE ST/END POINTS. WHERE IS 


c 



DEFINED BY SUBROUTINE DSTNCE. 


c 

XMAX 

: 

THE VALUE OF THE MAXIMUM DISTANCE AWAY FROM 

THE 

c 



ORIGIN ALONG AN OCTANT LINE. 


c 

XMIN 

: 

THE VALUE OF THE MINIMUM DISTANCE AWAY FROM 

THE 

c 



ORIGIN ALONG AN OCTANT LINE. 


c 

XNEW 

z 

THE X VALUE OF AN INTERSECTION POINT ON AN 


c 



OCTANT LINE. 


c 

XPTS 

z 

THE X VALUES IN CLOCKWISE ORDER OF THE POLYGON 

c 



BEING FORMED. 


c 

YNEW 

z 

THE Y INTERSECTION VALUE. (SEE XNEW) 


c 

YPTS 

z 

THE Y POINTS. (SEE XPTS) 



c 

c 

C COMMON VARIABLES: 

C NONE 

C 

C 

C FILES USED: 

C SCRATCH: 

C TAPE4 : COLLECTS SITE II CARDS FOR SUBSEQUENT 

C POST PROCESSING. 

C 
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SUBROUTINE DSTNCE 


C 
C 
C 

C OBJECTIVE: 

C TO COMPUTE THE DISTANCE FROM THE ORIGIN TO EVERY START 

C AND END POINT IN THE START/END POINT OCTANT TABLE AND 

C TO LOCATE TOE OCTANT ORDINAL FOR EACH START/END POINT. 

C 
C 

C METHOD: 

C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 


THE START/END POINT OCTANT TABLE. 

THE X POINTS OF THE CONTOUR BEING 
PROCESSED. 

THE Y POINTS OF THE CONTOUR BEING 
PROCESSED. 

THE NUMBER OF X,Y POINTS. 

THE TABLE OF START/END POINT DISTANCES 
TO THE ORIGIN. 

THE TABLE OF OCTANT ORDINALS FOR EACH 
START/END POINT. 

C 

C CALLING ROUTINE: 

C POLY 


A- 30 



GIVEN THE ORDINAL OF A START POINT WHICH IS INSIDE OR ON 
OCTANT LINE I, IT IS POSSIBLE TO CALCULATE A NEW POINT ON 
TOE OCTANT LINE BY FORMING THE INTERSECTION OF TOE OCTANT 
LINE AND THE LINE FORMED BY THE START POINT AND ITS PRED- 
ECESSOR. ONCE THE INTERSECTION POINT IS OBTAINED THE 
DISTANCE FORMULA IS USED WITH THE ORIGIN AS THE 2ND POINT 
TO FIND THE DISTANCE TO THE ORIGIN. THE INTERSECTION 
POINT IS ALSO PASSED TO FUNCTION IOCTNT WHICH RETURNS THE 
ORDINAL OF TOE OCTANT LINE THAT THE START/END POINT 
RESIDES ON. 


I. 



c 

c 

c 

c 

G 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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SUBPROGRAMS USED: 

INTSCT : 

CALCULATES TOE INTERSECTION OF A LINE WITH AN 
OCTANT LINE GIVEN 2 POINTS THAT ARE KNOWN TO 
STRADDLE AN OCTANT LINE. 

FUNCTIONS USED: 

IOCTNT : 

DETERMINES TOE ORDINAL OF THE OCTANT WHERE A 
POINT RESIDES. 

SQRT : 

FORTRAN SQUARE ROOT FUNCTION FOR REALS. 

LOCAL VARIABLE DICTIONARY: 

I : 

DO LOOP INDEX FOR OCTANT ORDINALS. 

INPT : 

A POINT WHICH IS LOCATED INSIDE A CERTAIN 
OCTANT. 

J : 

DO LOOP INDEX TO START/END POINTS. 

OUTPT : 

A POINT WHICH IS IN AN ADJACENT OCTANT RELATIVE 
TO EITHER ITS PREDECESSOR OR SUCCESSOR POINT. 

RAD : 

SPECIAL FORM OF THE DISTANCE FORMULA WHEN ONE 
POINT IS TOE ORIGIN: D = SQRT (RAD) WHERE 

RAD = X**2 + Y**2. 

XNEW : 

X COORDINATE OF INTERSECTION POINT WHOSE 
DISTANCE IS TO BE COMPUTED. 

YNEW : 

Y COORDINATE (SEE XNEW) . 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE INTSCT 
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c 

c 

c 

c 
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c 
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c 

c 

c 

c 

c 

c 

c 

c 
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OBJECTIVE: 

GIVEN TOO POINTS, CALCULATE THE EQUATION OF THE LINE THEY 
LIE ON, DETERMINE THE UNIQUE OCTANT OF EACH POINT TO FIND 
THE OCTANT LINE THEY STRADDLE AND FIND THE POINT OF 
INTERSECTION ON THE OCTANT LINE. 


METHOD: 

THE OCTANT LINES ARE THE SPECIAL FORMS: Y = X, Y = -X, 

Y = 0, X = 0. GIVEN THE 2 POINTS IT IS POSSIBLE TO 
DETERMINE THE STRADDLED OCTANT LINE USING THE FUNCTION 
IOCTNT WHICH DETERMINES THE OCTANT ORDINAL FOR A GIVEN 
POINT. ONCE THE OCTANT LINE IS KNOWN ITS EQUATION IS 
ALSO KNOWN. THE EQUATION OF THE LINE FOR THE 2 GIVEN 
POINTS IS OBTAINED USING THE POINT-SLOPE FORMULA. THEN 
CRAMER'S RULE IS USED TO CALCULATE THE INTERSECTION X, Y 
VALUES. SPECIAL CASES (INFINITE SLOPE) ARE HANDLED. 


PARAMETERS: 

INPUT: 

XIN : X VALUE OF INSIDE POINT. 

YIN : Y VALUE OF INSIDE POINT. 

XOUT : X VALUE OF OUTSIDE POINT. 

YOUT : Y VALUE OF OUTSIDE POINT. 

OUTPUT: 

X : X VALUE OF INTERSECTION POINT. 

Y : Y VALUE OF INTERSECTION POINT. 


CALLING ROUTINE: 
POLY 
DSTNCE 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 


C IOCTNT : DETERMINES ORDINAL OF OCTANT WHERE A PT RESIDES 

C 
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LOCAL VARIABLE DICTIONARY: 

A : THE X COEFFICIENT FOR THE EQUATION OF AN OCTANT 


B 

C 

D 

DELTA 

DX 

DY 

E 

F 

KI 

KJ 

LFLAG 


LSTOCT : 


M 


LINE. 

THE Y COEFFICIENT AS IN A ABOVE. 

THE CONSTANT AS IN A ABOVE. 

THE NEGATIVE SLOPE OF THE LINE FORMED BY 
(XIN, YIN) AND (XOUT, YOUT) . 

THE DENOMINATOR OR "DELTA" FOR CRAMER'S RULE. 
DETERMINENT AFTER SUBSTITUTING FOR X. 
DETERMINENT AFTER SUBSTITUTING FOR Y. 

Y COEFICIENT FOR THE LINE GIVEN 2 PTS. 

THE CONSTANT FOR THE LINE GIVEN 2 PTS. 

THE OCTANT ORDINAL FOR THE INSIDE POINT. 

THE OCTANT ORDINAL FOR THE OUTSIDE POINT. 
IDENTIFIES OCTANT LINE EQUATION WHEN SLOPE IS 
INFINITE FOR EQUATION OF A LINE. 

PARAMETER TO IOCTNT FUNCTION WHICH TELLS WHAT 
THE PREVIOUS OCTANT ORDINAL WAS. SET TO ZERO 
FOR THIS USAGE. 

THE SLOPE OF THE LINE FOR THE INSIDE/OUTSIDE 


POINTS. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE PLOT 


OBJECTIVE: 

TO INITIALIZE THE PLOT VECTOR FILE, CALCULATE THE SCALE 
FACTOR, DRAW THE COMPASS ROSE AND RESET THE PLOT ORIGIN. 


PARAMETERS: 

INPUT: 

XPTS 

YPTS 

N 

ITITLE 

LINE 

FACT 

OUTPUT: 

SCALE 


THE X VALUES OF THE POLYGON TO BE PLOTTED 
THE Y VALUES OF THE POLYGON TO BE PLOTTED 
THE NUMBER OF X, Y POINTS. 

THE 30 CHARACTER AIRPORT TITLE. 

THE LDN METRIC TITLE INFORMATION. 

HALF THE LENGTH OF THE X,Y AXIS. 

THE PLOTTING SCALE FACTOR. 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 

CALPLT : GRAPHICS ROUTINE FOR MOVING THE "PEN". 

PSEUDO : GRAPHICS ROUTINE FOR INITIALIZING THE PLOT 

VECTOR FILE (PVF) . 

ROSE : DRAWS THE COMPASS ROSE AND TITLES. 


FUNCTIONS USED: 

ABS : ABSOLUTE VALUE FOR REAL NUMBERS. 


LOCAL VARIABLE DICTIONARY: 

I : DO LOOP INDEX OVER THE RANGE OF X,Y POINTS. 

ZMAX : THE MAX VALUE OF ALL X,Y DATA POINTS. 

COMMON VARIABLES: 

NONE 

FILES USED: 

INPUT: 

TAPE 2 - SEQUENTIAL INPUT FILE USED WITH 

COUNTERMEASURES FOR READING IN THE 
PLOT SCALE OF THE BASELINE CASE 

OUTPUT: 

TAPE 2 - SEQUENTIAL FILE USED WITH THE BASELINE 
CASE FOR WRITING THE PLOT SCALE OF THE 
BASELINE CASE 
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SUBROUTINE ROSE 


C 
C 
C 

C OBJECTIVE: 

C TO PLOT THE COMPASS ROSE, PLOT BOUNDARIES AND LABELS. 

C 
C 

C METHOD: 

C 
C 
C 
C 
C 
C 

C PARAMETERS: 

C INPUT: 

C AXISLEN: 

C SCALE : 

C 

C ITITLE : 

C LINE : 

C 
C 

C CALLING ROUTINE: 

C PLOT 

C 


DRAWS BOUNDARY LINES (BOX) AND 3-LETTER 
COMPASS MNEUMONICS. 

GRAPHICS ROUTINE FOR PEN POSITIONING. 
PARTIALLY NOTATES THE PLOT LABELS. 
GRAPHICS CHARACTER PLOTTING ROUTINE. 


C FUNCTIONS USED: 
C NONE 

C 


.1 

it 


C SUBPROGRAMS USED: 
C BOUND : 
C 

C CALPLT : 
C LABELS : 
C NOTATE : 
C 


THE HALF-LENGTH OF THE X,Y AXES. 

THE SCALE FACTOR; USED TO DRAW THE SCALE 
CHART. 

THE 30-CHARACTER AIRPORT TITLE. 

THE LDN METRIC VALUE. 


DRAW A SQUARE USING THE GIVEN HALF-LENGTH FOR THE X, Y 
AXES; THEN DRAW + AND - X,Y AXES AND THE 2 LINES Y = X 
AND Y = -X (THE DIAGONALS OF THE SQUARE) . ANNOTATE EACH 
OCTANT WITH ITS 3-LETTER COMPASS MNEUMONIC. 


I 
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c 
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LOCAL VARIABLE DICTIONARY: 

CENTER : X COORDINATE FOR THE CENTER OF TOE SCALE CHART. 

HT : THE HEIGHT OF CHARACTERS IN FLOATING FT INCHES. 

IHOL : HOLLERITH SIRING ARRAY CONTAINING THE 3-LETTER 

COMPASS DIRECTION MNEUMONICS. 

S : THE DISTANCE FROM TOE LEFT MARGIN OF A SYMBOL 

TO THE LEFT MARGIN OF THE NEXT SYMBOL. 

TOTW : THE TOTAL WIDTH OF TOE SCALE LEDGER. 

W : THE WIDTH OF 1 CHARACTER OR SYMBOL. 

WADJ : THE WIDTH ADJUSTMENT WHICH IS USED TO 

COMPENSATE FOR CHARACTER HEIGHT. 

X : THE X COORDINATE PASSED TO THE BOUNDARY LINE 

ROUTINE. 

XLEN : THE LENGTH OF THE 2 MILE SCALE LINE IN INCHES. 

XI : THE LEFT MOST X COORDINATE OF THE SCALE LINE. 

X2 : THE RIGHT MOST COORDINATE OF TOE SCALE LINE. 

X3 : THE CENTER PT X COORDINATE OF THE SCALE LINE. 

Y : TOE COORDINATE PASSED TO THE BOUNDARY LINE 

ROUTINE. 

Y1 : THE Y COORDINATE OF TOE LEFT MOST PART OF TOE 

SCALE LINE. 

Y2 : TOE COORDINATE OF TOE POSITIVE END OF TOE 

VERTICAL SCALE LINE TIC MARKS. 

Y3 : LIKE Y2 ONLY THE NEGATIVE END. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE LABELS 


OBJECTIVE: 

TO DRAW THE AIRPORT NAME AND "LDN FOOTPRINT" AT THE TOP 
OF THE PLOT; TO ESTABLISH VALUES, LATER USED BY RESTPLT, 
TO COMPLETE THE PLOT LABELING. 


PARAMETERS: 

INPUT: 

CENTER : 


ITITLE 

LDNST 

CALLING ROUTINE: 

ROSE 


THE CENTER X AND Y VALUE OF THE COMPASS 
ROSE. 

THE AIRPORT NAME 

THE STARTING LDN METRIC. 


SUBPROGRAMS USED: 

NOTATE : GRAPHICS ROUTINE USED TO DRAW CHARACTER STRINGS 


FUNCTIONS USED: 

FLOAT : INTEGER TO REAL CONVERSION 
OR : LOGICAL OR OF 2 OR MORE ARGUMENTS. 

SHIFT : LEFT CIRCULAR ( LOGICAL) /RT ARITHMETIC WITH SIGN 

EXTENSION SHIFT. 


LOCAL VARIABLE DICTIONARY: 

HT : THE HEIGHT IN INCHES OF CHARACTERS TO BE 

NOTATED. 

I : DO LOOP INDEX USED FOR FINDING A.P. NAME WIDTH. 

ICHAR : COMPUTED CHARACTER PACKED INTO THE RANGE OF 

LDN CONTOUR TITLE. 

IDUM : DUMMY LIST VARIABLE USED FOR DECODE. 

IREM : VARIABLE USED TO DO DIGIT BY DIGIT ENCODE OF 

LDN METRIC VALUES FOR LDN RANGE TITLE. 

ISTRIN : CHARACTER STRING (2 WORDS) WITH HOLES FOR LOW 

AND HIGH LEW RANGE VALUES: 

" 55 DB TO 85 DB " IS THE IDEAL FINAL FORM. 

J : BACKWARD COUNTER FROM 80 AND 40 IN CONTRAST WITH I. 

LINE : DECODE TARGET (8 WORDS ONTO 80) FOR A.P. TITLE. 

LINE2 : 40-WORD ARRAY CONTAINING CASE TITLE 
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MET ' : THE ALPHA FORM OF THE LDN METRIC VALUE IN F5.1 

METST : THE INTEGER FORM OF MET. (THE STARTING VALUE) . 

S : THE DISTANCE FROM THE LEADING EDGE OF A SYMBOL 

TO THE LEADING EDGE OF THE NEXT SYMBOL. 

TOTW : THE TOTAL WIDTH OF A CHARACTER STRING IN INCHES 

W : THE WIDTH OF ONE SYMBOL. 

WADJ : WIDTH ADJUSTMENT FACTOR BASED ON SYMBOL HEIGHT. 

X : THE X COORDINATE OF THE LOWER LEFT HAND CORNER 

OF THE 1ST CHARACTER OF A STRING TO BE DRAWN 
SUCH THAT IT IS CENTERED. 

XMET : THE REAL FORM OF THE LDN STARTING CONTOUR. 

Y : Y COORDINATE - SEE X. 


COMMON VARIABLES : 

/ITITL/ 

ICASE1 : 4-WORD ARRAY CONTAINING CASE TITLE 

/PLTINFO/ 

CENT : THE CENTER X AND Y VALUE OF THE COMPASS 
ROSE 

HT : THE HEIGHT IN INCHES OF CHARACTERS TO BE 
NOTATED 

ISTRIN : CHARACTER STRING (2 WORDS) USED TO NOTATE 
RANGE OF LDN VALUES 

S : THE DISTANCE FROM THE LEADING EDGE OF A 

SYMBOL TO THE LEADING EDGE OF THE NEXT 
SYMBOL 

W : REAL VARIABLE REPRESENTING THE WIDTH OF 

ONE SYMBOL 


FILES USED: 

NONE 


A- 3 8 



SUBROUTINE BOUND 
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OBJECTIVE: 

TO DRAW THE PLOTTING BOUNDARY LINES AND ANNOTATE THE 
COMPASS ROSE DIRECTION MNEUMONICS. 


PARAMETERS: 

INPUT: 

AXLEN : HALF TOE LENGTH OF TOE X,Y AXES. 

X : THE X PT OF ORIGIN FOR THE BOUNDARY LINE 

TO BE DRAWN. 

Y : Y COORDINATE (SEE X) . 

INDEX : THE INDEX TO THE 3-LETTER DIRECTION 

MNEUMONICS TO BE NOTATED ON THIS CALL. 
IANGL : 0 FOR VERTICAL BOUNDARY, 1 FOR HORIZONTAL. 

IHOL : THE ARRAY OF 3 LETTER COMPASS HEADING 

MNEUMONICS. 


CALLING ROUTINE: 
ROSE 


SUBPROGRAMS USED: 

CALPLT : X,Y PEN MOVEMENT FOR GRAPHICS. 

NOTATE : NOTATES CHARACTER STRINGS ON GRAPHIC DISPLAY. 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

HT : CHARACTER HEIGHT. 

THETA : ANGLE OF ORIENTATION FOR CHARACTER STRINGS. 

TOIWDT : TOTAL WIDTH OF 3 LETTER MNEUMONIC. 

WADJ : CHARACTER WIDTH ADJUSTMENT FACTOR. 

X4 : THE DESTINATION PT FOR PLOT PEN IN DOWN POS. 

X COORD WHEN IANGL = 1, Y COORD WHEN IANGL = 0. 
X6 : X COORD OF LOWER LEFT HAND CORNER OF 1ST 

MNEUMONIC PLOTTED. 

X7 : SAME AS X6 BUT FOR 2ND MNEUMONIC PLOTTED (THERE 

ARE 2 MNEUMONICS FOR EACH BOUNDARY LINE. 

X8 : THE X COORDINATE FOR THE HORIZONTAL BOUNDARY 

LINE MNEUMONICS. 

Y3 : THE Y COORDINATE OF HORIZONTAL BOUNDARY LINE 

MNEUMONICS. 

Z1 : THE X,Y QUADRANT MIDPOINT (ONE FOURTH THE 

LENGTH OF TOE X, Y AXES. 
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OBJECTIVE: 

TO DRAW A POLYGON GIVEN THE X, Y POINTS AND A SCALE FACTOR 


METHOD: 

CONVERT THE SCALE FACTOR TO MILES UNITS AND DIVIDE ALL 
DATA VALUES BY NEW SCALE FACTOR. THEN USE MULTIPLE POINT 
DRAW GRAPHICS ROUTINE. 


PARAMETERS: 

X 

Y 

N 

SCALE 


X POINTS FOR POLYGON TO BE PLOTTED. 
Y POINTS (SEE X) . 

NUMBER OF X, Y POINTS IN POLYGON. 

THE SCALE FACTOR. 


CALLING ROUTINE: 
POLY 


SUBPROGRAMS USED: 
DRAW : 


DRAWS A LINE GIVEN A SERIES OF SCALED DATA PTS. 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE 
DIV 
I 

XS 

YS 


DICTIONARY: 

: SCALE FACTOR AFTER CONVERSION TO STATUTE MILES. 

: DO LOOP VARIABLE OVER RANGE OF POINTS. 

: SCALED X POINTS. 

: SCALED Y POINTS. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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SUBROUTINE BLDARR 
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OBJECTIVE: 

TO EVALUATE TOE SITE II INPUT DECK, SUBDIVIDING ANY DECK 
WITH MORE THAN 150 POINTS INTO 2 SITE II RUNS; TO PUT 
END-OF-FILE MARKS AT TOE END OF INPUT DECKS FOR EACH 
SITE II RUN. ADDITIONALLY, TO COMPUTE VALUES FOR REGISTER 
SETTINGS, WITH NCON INDICATING THE TOTAL NUMBER OF 
SITE II RUNS, AND NSTCON PROVIDING AN INDEX TO THE LDN 
VALUE OF TOE FIRST CONTOUR FOUND. 


METHOD: 

FOR EACH SITE II RUN (INDICATED BY AN 'A' CARD), COUNT 
THE CORRESPONDING NUMBER OF POINTS ON THE ASSOCIATED 
'D' CARDS; STORE TOE TOTAL NUMBER OF POINTS IN ARRAY 
IPTS; FOR EACH RUN EXCEEDING 150 POINTS, THE ASSOCIATED 
GROUPS OF 'C* AND 'D' CARDS ARE SUBDIVIDED. 


PARAMETERS: 

INPUT: 

ICASE: 8-WORD ARRAY WITH WORDS 5 THROUGH 8 
CONTAINING THE CASE NAME 

OUTPUT: 

NCON: INTEGER VARIABLE, TO BE USED FOR REGISTER 
SETTING, INDICATING THE TOTAL NUMBER OF 
SITE II EXECUTIONS REQUIRED 
NSTCON: INTEGER VARIABLE, TO BE USED FOR REGISTER 
SETTING, THAT IS AN INDEX TO THE LDN VALUE 
OF THE FIRST CONTOUR FOUND 


CALLING ROUTINE: 

AUPHS1 (MAIN PROGRAM) 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 
NONE 
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LOCAL VARIABLE DICTIONARY: 

IA: INTEGER ARRAY USED TO CONTAIN THE CARD IMAGE OF A 
SITE II 'A' CARD 

IARAY: 2-DIMENSION INTEGER ARRAY, WITH THE FIRST DIMEN- 
SION USED TO REPRESENT A SINGLE LDN CONTOUR VALUE, 
AND THE SECOND DIMENSION USED TO INDICATE THE NUM- 
BER OF POINTS IN A SINGLE POLYGON FOR THE ASSOCIATED 
LDN CONTOUR VALUE 

IB: INTEGER ARRAY USED TO CONTAIN THE CARD IMAGE OF A 
SITE II 'B' CARD 

IBLK: INTEGER ARRAY USED TO CONTAIN THE FIRST SITE II 
CARD IMAGE 

ICHAR: INTEGER VARIABLE CONTAINING, IN THE LAST 6 BITS, 
THE FIRST 6 BITS OF A SITE II CARD IMAGE 
ICHAR1: INTEGER VARIABLE CONTAINING, IN THE LAST 6 BITS, 
THE FIRST 6 BITS OF A SITE II CARD IMAGE 
ICHAR2: INTEGER VARIABLE CONTAINING, IN THE LAST 6 BITS, 
THE FIRST 6 BITS OF A SITE II CARD IMAGE 
ICOUNT: INTEGER VARIABLE USED TO INDICATE A SPECIFIC 
POLYGON ASSOCIATED WITH A SPECIFIC LDN CONTOUR 
ICTR: INTEGER VARIABLE USED TO INDICATE THE NUMBER OF 
POINTS IN A PARTICULAR POLYGON 
IF ACT: INTEGER VARIABLE INDICATING THE TOTAL NUMBER OF 
LDN CONTOURS FOUND 

IJ: INTEGER VARIABLE IB ED AS POLYGON INDEX 
IJK: INTEGER VARIABLE USED AS LDN CONTOUR INDEX 
IPTS: INTEGER ARRAY USED TO CONTAIN THE TOTAL NUMBER OF 
POINTS FOR EACH LEW CONTOUR 
11: INTEGER VARIABLE, CONTAINING IN THE LAST 6 BITS, 

THE FIRST 6 BITS OF A SITE II CARD IMAGE 
ISA: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I5B: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I6A: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I6B: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I7A: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I8A: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

I8B: INTEGER USED AS TEMP. STORAGE FOR PART OF CASE TITLE 

LINE: INTEGER ARRAY USED TO CONTAIN A SITE II CARD 
IMAGE 

SPACE: INTEGER ARRAY USED TO CONTAIN THE FIRST FOUR 
WORDS OF A SITE II CARD IMAGE 
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c 

C COMMON VARIABLES: 

C NONE 

C 

C 

C FILES USED: 

C INPUT: 

C NONE 

C OUTPUT: 

C TAPE1 - SEQUENTIAL FILE CONTAINING SITE II INPUT 

C DECK, EACH DESCRIBING 150 POINTS OR LESS, 

C WITH INPUT DECKS SEPARATED BY END-OF-FILE MARKS 

C SCRATCH: 

C TAPE4 - SEQUENTIAL FILE CONTAINING CARD IMAGES OF 

C UNPROCESSED SITE II INPUT DECKS 

C 
C 
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SUBROUTINE RESTPLT 
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OBJECTIVE: 

SUBROUTINE RESTPLT PLOTS A LABEL INDICATING THE RANGE 
OF DB VALUES WHOSE CONTOUR POINTS ARE PLOTTED; ADDITION- 
ALLY, IF MORE THAN ONE CONTOUR IS PLOTTED, THE STEP SIZE 
IS NOTATED ON THE PLOT 


PARAMETERS: 

INPUT: 

FACT : TOE LENGTH OF TOE POSITIVE X AND Y AXES IN 
INCHES 

MET : NOISE METRIC IN CHARACTER FORM 
RUNCNT: INTEGER VARIABLE USED TO COUNT CONTOURS 
OUTPUT: 

NONE 


CALLING ROUTINE: 
AUPHS1 


SUBPROGRAMS USED: 

CALPLT : LRCGOSF ROUTINE USED TO RESET ORIGIN 
NOTATE : LRCGOSF ROUTINE USED TO PLOT A CHARACTER 
STRING 


FUNCTIONS USED: 

AND : STANDARD FTN AND-ING FUNCTION 
OR : STANDARD FTN OR-ING FUNCTION 
SHIFT : STANDARD FTN SHIFT-ING FUNCTION 


LOCAL VARIABLE DICTIONARY: 

CENTER : REAL VARIABLE REPRESENTING THE CENTER X AND Y 
VALUE OF THE COMPASS ROSE 

ICHAR : VARIABLE USED IN THE FORMATION OF THE UPPERMOST 
LDN CONTOUR VALUE PLOTTED 

IREM : VARIABLE USED IN THE FORMATION OF TOE UPPERMOST 
LDN CONTOUR VALUE PLOTTED 

ITEMP1 : VARIABLE USED IN THE FORMATION OF CHARACTER 
STRING WHEN ONLY ONE CONTOUR IS PLOTTED 

ITEMP2 : VARIABLE USED TO STORE THE INTEGER VALVE OF 
THE UPPERMOST LDN VALUE 

TOTW : REAL VARIABLE REPRESENTING THE LENGTH IN INCHES 
OF THE CHARACTER STRING TO BE PLOTTED 
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X : REAL VARIABLE REPRESENTING THE X-COORDINATE OF 

THE STARTING POINT OF THE CHARACTER STRING TO 
BE NOTATED 

Y : REAL VARIABLE REPRESENTING THE Y-COORDINATE OF 

THE STARTING POINT OF THE CHARACTER STRING TO 
BE NOTATED 


COMMON VARIABLES: 

/PLTINFO/ 

CENT : THE CENTER X AND Y VALUE OF THE COMPASS 
ROSE 

HT : THE HEIGHT IN INCHES OF CHARACTERS TO BE 
NOTATED 

ISTRIN : CHARACTER STRING (2 WORDS) USED TO NOTATE 
RANGE OF LDN VALUES 

S : THE DISTANCE FROM THE LEADING EDGE OF A 

SYMBOL TO THE LEADING EDGE OF THE NEXT 
SYMBOL 

W : REAL VARIABLE REPRESENTING THE WIDTH OF 

ONE SYMBOL 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

NONE 

SCRATCH: 

NONE 
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OBJECTIVE: 

PROGRAM AUREPT CONSOLIDATES THE DEMOGRAPHIC VARIABLES 
FROM A SERIES OF SITE II DEMOGRAPHIC PROFILE REPORTS FOR 
AN AIRPORT COMMUNITY WITHIN THE 65 DB LDN NOISE CONTOUR. 
TEN PRINTED REPORTS ARE GENERATED; ONE REPORT FOR EACH 
OCTANT OF THE IMPACTED AREA, A SUMMARY OF ALL OCTANT 
REPORTS AND A NOISE IMPACT SUMMARY FOR THE IMPACTED AREA. 


METHOD: 

INPUT TO PROGRAM AUREPT CONSISTS OF A SITE II 
DEMOGRAPHIC PROFILE REPORT FOR EACH OF A SERIES OF 
POLYGONAL AREAS (PA) . EACH PA REPRESENTS THE AREA 
ENCLOSED BY THE COORDINATES OF A CONSTANr LEVEL NOISE 
CONTOUR INTERSECTED WITH TWO COMPASS ROSE OCTANT LINES. 
THE DEMOGRAPHIC VARIABLES FROM THE PROFILE REPORTS ARE 
SUMMED FOR EACH OCTANT OF EACH CONTOUR. DIFFERENCES ARE 
CALCULATED BETWEEN EQUIVALENT OCTANTS OF ADJACENT 
CONTOURS SUCH THAT DEMOGRAPHIC VARIABLES ARE STORED FOR 
THE AREAS (BANDS) BETWEEN EACH CONTOUR FOR EACH OCTANT. 

THE DEMOGRAPHIC VARIABLES ARE ACCUMULATED IN A 
3-DIMENSIONAL ARRAY SUCH THAT EACH PLANE REPRESENTS AN 
OCTANT WITH THE 9TH PLANE REPRESENTING THE SUM OF ALL 
OCTANTS. THE ROWS OF EACH PLANE REPRESENT THE 
DEMOGRAPHIC VARIABLES WHILE THE COLUMNS REPRESENT THE 
CONTOUR BANDS. THE AGGREGATED POPULATION DATA AS SUMMED 
AND DIFFERENCED FROM THE DEMOGRAPHIC DATA IS USED TO 
CALCULATE A NOISE IMPACT SUMMARY WHICH INCLUDES THE 
CALCULATION OF A LEVEL WEIGHTED POPULATION, NOISE IMPACT 
INDEX, EQUIVALENT NOISE LEVEL, PERCENT HIGHLY ANNOYED AND 
POPULATION HIGHLY ANNOYED. 


SUBPROGRAMS USED: 


BLDTAB : 

READS THE PA REPORTS, CALCULATES THE BAND DATA 
AND STORES IT IN THE 3-DIMENSIONAL TABLE. 

CHART : 

OUTPUTS THE OCTANT REPORTS ON PRINTER PAPER 
USING THE CONTENTS OF THE 3-DIMENSIONAL TABLE. 

HEADER : 

PERFORMS PAGE EJECT AND PRINTS PAGE HEADINGS 
AND TITLES FOR THE OCTANT REPORTS. 

SUMARY : 

CALCULATES AND PRINTS TOE VARIABLES FOR TOE 
NOISE IMPACT SUMMARY. 

FUNCTIONS USED: 


NONE 
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LOCAL VARIABLE DICTIONARY : 

ALINE : 2-DIMENSIONAL ARRAY USED TO STORE AREAS 

OF 65, 75, 85 DB CONTOURS. 

GR2 : USER-DESIGNATED GROWTH RATE 

I : USED AS A DO LOOP INDEX IN 50 LOOP AND 100 LOOP 

IN THE 50 LOOP IT INDEXES THE ARRAY TABLE FOR 
PRESET TO ZEROES AND IN THE 100 LOOP IT IS THE 
OCTANT COUNTER. 

ICASE : 60-WORD ARRAY CONTAINING CASE TITLE 

IR3 : VALUE OF REGISTER 3, AS DETERMINED BY GETCON 

ISIG1 : FLAG USED TO INDICATE THE ABSENCE OF THE 75 DB 

(OR GREATER) CONTOUR 

ISIG2 : FLAG USED TO INDICATE THE ABSENCE OF THE 65 DB 

(OR GREATER) CONTOUR 

IYEAR : USER-DESIGNATED YEAR FOR POPULATION UPDATES 

J : USED AS A DO LOOP INDEX IN 75 LOOP. IT IS THE 

COLUMN (BAND) POINTER IN THE ARRAY TABLE. 

KADRS : AN 80 WORD ARRAY THAT STORES THE AIRPORT 

ADDRESS IN 80R1 FORMAT. 

KEMPTY : A COUNTER WHOSE VALUE DETECTS THE PRESENCE OF 

AN OCTANT WITH NO POPULATION IN ANY BAND. 

LDNS : A 8 WORD ARRAY THAT STORES THE ASCENDING LDN 

VALUES FOR EACH CONTOUR. DEFINED IN SUBROUTINE 
BLDTAB. 

MINUS : INTEGER VARIABLE USED FOR LANDUSE COUNTER- 

MEASURES TO DETERMINE DB INSULATION FACTOR 

NAME : AN 80 WORD ARRAY THAT STORES THE AIRPORT NAME 

IN 80R1 FORMAT. 

NCON : THE INDEX TO THE STARTING LDN CONTOUR. DEFINED 

IN SUBROUTINE BLDTAB BY SAMPLING THE JOB 
CONTROL REGISTER, R2. 

TABLE : THE 3-DIMENSIONAL TABLE BUILT UP BY SUBROUTINE 

BLDTAB. 

COMMON BLOCK CROSS-REFERENCE LIST: 

NONE 


COMMON VARIABLES: 
NONE 
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FILES USED: 

INPUT: 

TAPE 5 - SEQUENTIAL INPUT FILE CONSISTING OF ALL THE 
SITE II PRINTED OUTPUT. FORMATTED MODE. 

TAPE1 - SEQUENTIAL INPUT FILE CONSISTING OF AIRPORT 
DESCRIPTIVE DATA. FORMATTED MODE. 

CARDS 1-6: INFORMATION USED BY ALAMO PROGRAMS THAT 
SHARE THIS FILE 

CARD 7: LATITUDE, LONGITUDE, DISPLACEMENTS. 

COL1 - COL2 = LATITUDE DEGREES 
COL4 - COL5 = LATITUDE MINUTES 
COL7 - COL8 = LATITUDE SECONDS 
COL11 - COL13 = LONGITUDE DEGREES 
COL15 - COL16 = LONGITUDE MINUTES 
COL 18 - COL19 = LONGITUDE SECONDS 
COL25 = N OR S (NORTR/SOUTH) 

COL 2 6 - COL30 = F5.2 WITH NO SIGN 
COL32 = E OR W (EAST/WEST) 

COL33 - COL37 = F5.2 WITH NO SIGN 

CARD 8: 2F10.2 - X AND Y AXES TRANSLATION 
CONSTANTS. 

CARD 9: COL1 - COL5: F5.2 - SCALE FACTOR FOR CONTOUR 
GRAPHIC EQUALS HALF TOE LENGTH OF X AND Y AXES 
COL7 - COLIO: YEAR FOR POPULATION UPDATES 
COL11 - COL16 : USER DESIGNATED GROWTH RATE 
EQUALS HALF THE LENGTH OF X AND Y AXES. 

CARD 10: AIRPORT NAME IN 8A10 FORMAT. 

CARD 11: AIRPORT ADDRESS IN 8A10 FORMAT. 

OUTPUT: 

TAPE6 - LINE EH INTER FILE FOR ALL 10 REPORTS. 

TAPE2 - LINE PRINTER FILE FOR NOISE IMPACT SUMMARY 
REPORT ONLY 

SCRATCH: 

NONE 
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OBJECTIVE: 

TO FIND AND READ THE NEXT SITE II POLYGON REPORT ON TAPE 5 


METHOD: 

RDREPT READS TAPE5, LINE BY LINE, UNTIL IT FINDS THE NEXT 
OCCURRENCE OF THE CHARACTER SIRING " PROFILE R" STARTING 
IN COLUMN 41. WHEN THE STRING IS FOUND RDREPT WILL BE 
POSITIONED AT THE TOP OF A SITE II DEMOGRAPHIC PROFILE 
REPORT. THE DEMOGRAPHIC VARIABLES FOR THE NEXT POLYGON 
ARE THEN READ THROUGH A SERIES OF FORMATTED READS. 

ASSUMPTIONS: 

RDREPT ASSUMES THAT THE SITE II REPORTS ARE IN THE PROPER 
ORDER FOR LATER PROCESSING. IT MAKES NO ATTEMPT TO 
RECOGNIZE EITHER THE OCTANT NUMBER OR THE LDN VALUE. 

PARAMETERS: 

INPUT: 

GR2 : USER-DESIGNATED GROWTH RATE 

IYEAR : USER-DESIGNATED YEAR OF POPULATION UPDATE 

OUTPUT: 

IR3 : VALUE OF REG. 3 AS DETERMINED BY GETCON 

ISIG1 : FLAG USED TO INDICATE ABSENCE OF 75 E 8 
(OR GREATER) CONTOUR 

ISIG2 : FLAG USED TO INDICATE ABSENCE OF 65 IB 

(OR GREATER) CONTOUR 

CALLING ROUTINE: 

BLDTAB 

SUBPROGRAMS USED: 

GETCON : COMPASS ROUTINE THAT SAMPLES REGISTERS 


FUNCTIONS USED: 

EOF: END OF FILE CHECK ON TAPE5. 


LOCAL VARIABLE DICTIONARY: 

I : USED AS INDEX VARIABLE IN IMPLIED DO I/O LISTS. 

IR1 : INTEGER VARIABLE REPRESENTING VALUE OF REGISTER 1 

IR2 : INTEGER VARIABLE REPRESENTING VALUE OF REGISTER 2 

IR3 : INTEGER VARIABLE REPRESENTING VALUE OF REGISTER 3 

KEY : STORES 10 CHARACTER STRING WHICH UNIQUELY 

IDENTIFIES A SITE II DEMOGRAPHIC PROFILE REPORT 
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c 

COMMON VARIABLES: 




c 

/COMIV 




c 

ACHH 


THE NUMBER OF HOUSEHOLDS WITH 

c 



AIR CONDITIONING. 

c 

AREA 

2 

THE AREA IN SQUARE MILES OF THE POLYGON 

c 



BEING 

READ. 

c 

COL16P 

: 

NOT USED 

c 

EDUC 

• 

• 

5 WORD VECTOR STORING THE EDUCATED ADULTS 

c 



GREATER THAN 25 COLUMN INFO AS FOLLOWS: 

c 



EL 1 

0-8 GRADES 

c 



EL 2 

9-11 GRADES 

c 



EL 3 

12 TH GRADE 

c 



EL 4 

13-15 GRADES 

c 



EL 5 

16+ GRADES 

c 

FAMPOP 

• 

THE NUMBER OF PEOPLE WITH FAMILY STATUS. 

c 

FEM 

: 

6 WORD VECTOR STORING THE NUMBER OF WOMEN 

c 



IN THE FOLLOWING AGE CATEGORIES: 

c 



EL 1 

18 - 20 

c 



EL 2 

21 - 29 

c 



EL 3 

30 - 39 

c 



EL 4 

40 - 49 

c 



EL 5 

50 - 64 

c 



EL 6 

65+ 

c 

EMIN 


7 WORD VECTOR STORING THE NUMBER OF PEOPLE 

c 



IN THE FOLLOWING FAMILY INCOME BRACKETS: 

c 



EL 1 

$ 0 - 5 (THOUSANDS) 

c 



EL 2 

$ 5 - 7 

c 



EL 3 

$ 7 - 10 

c 



EL4 

$10 - 15 

c 



EL5 

$15 - 25 

c 



EL6 

$25 - 50 

c 



EL7 

$50+ 

c 

HOMVAL 

• 

7 WORD VECTOR STORING THE NUMBER OF PEOPLE 

c 



IN THE FOLLOWING HOME VALUE BRACKETS: 

c 



ELI 

$ 0 - 10 (THOUSANDS) 

c 



EL2 

$10 - 15 

c 



EL3 

$15 - 20 

c 



EL4 

$20 - 25 

c 



EL5 

$25 - 35 

c 



EL6 

$35 - 50 

c 



EL7 

$50+ 

c 

I OCT 

• 

• 

THE OCTANT NUMBER OF THE POLYGON BEING 

c 



READ. 


c 

KEOF 

• 

END OF FILE FLAG - SET TO 1 TO INDICATE 

c 



AN EOF WAS ENCOUNTERED. SENSED BY THE 

c 



CALLING ROUTINE. 

c 

LINE 

i 

AN 8 

WORD VECTOR FOR STORING 80 CHARACTER 

c 



LINE 

IMAGES OF THE SITE II REPORTS. 

c 

MALE 

: 

6 WORD VECTOR STORING MALE AGE INFO AS 

c 



DESCRIBED UNDER VARIABLE FEM. 
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1 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MGRPRO : 


NOHH : 
POP70 : 
POP81 : 
SNGEWL : 

TOTCWN : 
TOTRNT : 
TVHH : 
XLDN : 

ZERO : 


9 WORD VECTOR STORING THE NUMBER OF PEOPLE 

IN THE FOLLOWING OCCUPATIONAL CATEGORIES: 

ELI MANAGERIAI/PROFESSIONAL 

EL2 SALES 

EL3 CLERICAL 

EL4 CRAFT 

EL5 OPERATIVES 

EL6 LABORER 

EL7 FARM 

EL8 SERVICE 

EL9 PRIVATE 

THE NUMBER OF HOUSEHOLDS. 

THE 1970 CENSUS POPULATION FOR THIS POLY. 
THE 1981 CENSUS POPULATION FOR THIS POLY. 
THE NUMBER OF PEOPLE IN SINGLE FAMILY 
DWELLINGS. 

THE TOTAL NUMBER OF HOMEOWNERS. 

NOT USED. 

THE NUMBER OF HOUSEHOLDS WITH TELEVISION. 
THE NOISE CONTOUR METRIC (DBLDN) FOR THIS 
POLYGON. 

A FLAG SET TO 1 TO INDICATE A ZERO 
POPULATION POLYGON. 


FILES USED: 

INPUT: 

TAPE 5 - SEQUENTIAL INPUT FILE CONTAINING ALL OF THE 
SITE II OUTPUT FOR AN ALAMO CASE RUN. THIS 
CONSISTS OF A REPITITION OF THE FOLLOWING: 

1. SITE II REFLECTION OF A,B,C AND D CARDS 
FOR A PARTICULAR RUN. 

2. REFLECTION OF THE SORT/MERGE DIRECTIVES 
USED BY SITE II. 

3. SITE II DEMOGRAPHIC PROFILE REPORTS FOR 
INPUT DESCRIBED IN 1 ABOVE. 
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SUBROUTINE BLDTAB 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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OBJECTIVE: 

TO BUILD A 3-DIMENSIONAL TABLE WHICH STORES THE VALUES OF 
THE 16 DEMOGRAPHIC VARIABLES TO BE REPORTED FOR EACH 
OCTANT OF EACH LDN BAND. 


METHOD: 

BLDTAB CALLS SUBROUTINE RDREPT TO READ EACH SUCCESSIVE 
POLYGON REPORT. RDREPT DEFINES AND REDEFINES THE 
DEMOGRAPHIC VARIABLES IN NAMED COMMON BLOCK, COMM. 

BLDTAB ACCUMULATES THE VALUES FOR ALL POLYGONS OF EACH 
OCTANT FOR EACH LDN VALUE AND STORES THE ACCUMULATIONS IN 
LOCAL VARIABLES. WHEN ALL POLYGONS HAVE BEEN READ AND 
VALUES ARE LOCALLY STORED, BLDTAB BEGINS TO BUILD THE 
3-DIMENSIONAL TABLE BY FINDING THE DIFFERENCES OF THE 
ACCUMULATED VALUES OF THE DEMOGRAPHIC VARIABLES BETWEEN 
LDN CONTOURS. THIS IS DONE FOR EACH OCTANT. THE ROW 
DIMENSION OF THE TABLE CORRESPONDS TO THE DEMOGRAPHIC 
VARIABLES; THE COLUMN DIMENSION OF THE TABLE CORRESPONDS 
TO THE LDN CONTOUR BANDS AND THE PLANE DIMENSION 
CORRESPONDS TO THE 8 COMPASS ROSE OCTANTS. THE 9TH PLANE 
REPRESENTS A SUMMARY OF ALL THE OCTANTS. THE CAPACITY OF 
THE TABLE ACCOUNTS FOR THE FOLLOWING BANDS: 
65-70,70-75,75-80, 80-85 AND OVER 85. 

IT IS POSSIBLE THAT LOW END AND/OR TOP END LEW CONTOURS 
WILL BE MISSING. IN ANY CASE, THE BANDS WILL BE STORED 
IN THE UPPER DIMENSION OF THE COLUMNS OF THE TABLE. THIS 
MEANS THAT COLUMN 5 WILL ALWAYS BE THE "OVER XX" BAND 
WHERE XX CORRESPONDS TO THE HIGHEST LEW CONTOUR FOUND BY 
THE INTEGRATED NOISE MODEL. BLDTAB SAMPLES TOE JOB 
CONTROL REGISTER, R2 FOR THE STARTING COLUMN INDEX TO 
ACCOUNT FOR MISSING CONTOURS. R2 IS SET BY A SEPARATE 
PROGRAM WHICH DEFINES THE POLYGON AREAS AND THEREFORE 
KNOWS THE LDN RANGE OF VALUES. 

RESTRICTIONS AND ASSUMPTIONS: 

BLDTAB ASSUMES THAT: 

1. REGISTER R2 HAS BEEN SET PREVIOUS TO PROGRAM 
STARTUP. 

2. THE POLYGON REPORTS ARE IN ASCENDING ORDER OF 
LDN VALUE. 
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PARAMETERS: 

INPUT: 

GR2 : 

IYEAR : 
OUTPUT: 

ALINE : 

IR3 : 

ISIG1 : 

ISIG2 : 

TABLE : 


LDNS : 


ICON : 


USER-DESIGNATED GROWTH RATE 
USER-DESIGNATED YEAR OF POPULATION UPDATE 

2-DIMENSIONAL ARRAY CONTAINING AREAS OF 
65, 75, 85 EB CONTOURS. 

VALUE OF REG. 3 AS DETERMINED BY GETCON 
FLAG USED TO INDICATE ABSENCE OF 75 DB 
(OR GREATER) CONTOUR 
FLAG USED TO INDICATE ABSENCE OF 65 DB 
(OR GREATER) CONTOUR 

THE 3 DIMENSIONAL TABLE - TABLE (16, 5, 9) 
EACH PLANE REPRESENTS AN OCTANT EXCEPT THE 
9TH PLANE WHICH REPRESENTS THE SUMMARY OF 
ALL THE OCTANTS. THE 5 COLUMNS OF EACH 
PLANE CORRESPOND TO THE CONTOUR BANDS AS 
FOLLOWS: 

65-70 70-75 75-80 80-85 85+ 

EACH ROW CORRESPONDS TO A DEMOGRAPHIC 
VARIABLE AS FOLLOWS: 

1. SELF-NOISE (AMBIENT NOISE IN EBLDN) 

2. IMPACTED POPULATION 

3. AREA IN SQUARE KILOMETERS 

4. TOTAL POPULATION PER SQUARE KILOMETER 

5. AVERAGE GROWTH RATE, APR 

6. PERCENT FAMILY POPULATION 

7. AVERAGE AGE, ADULTS OVER 17 

8. PERCENT AGE OVER 65 

9. PERCENT OF PEOPLE WITH 16+ YEARS EDUC. 

10. PERCENT MANAGER/PROFESSIONALS 

11. AVERAGE FAMILY INCOME 

12. PERCENT SINGLE FAMILY DWELLINGS 

13. PERCENT HOME OWNERS 

14. AVERAGE HOME VALUE 

15. PERCENT HOUSEHOLDS WITH AIRCONDITIONER 

16. PERCENT HOUSEHOLDS WITH TV 

A 6 WORD ARRAY WHICH STORES THE 5 
POSSIBLE LDN CONTOUR VALUES INTENDED TO BE 
65 THROUGH 85 IN 5DB INTERVALS. THE 6TH 
WORD STORES THE MAXIMUM LEW CONTOUR FOUND 
WHICH WILL BE THE CONTENTS OF WORD 5. 
STORAGE FOR THE VALUE OF JOB CONTROL 
REGISTER, R2. DEFINED BY COMPASS 
SUBROUTINE GETCON. 
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CALLING ROUTINE: 
AUREPT 


SUBPROGRAMS USED: 

GETCON : READS JOB CONTROL REGISTERS. 

RDREPT : READS THE NEXT POLYGON REPORT ON THE INPUT FILE 


FUNCTIONS USED: 

ALQG10 : CALCULATES THE LOGARITHM, BASE 10 


LOCAL VARIABLE DICTIONARY: 

AC : AN 8X6 MATRIX WHOSE ROWS CORRESPOND TO THE 8 

OCTANTS AND WHOSE COLUMNS CORRESPOND TO THE 5 
POSSIBLE LDN VALUES. COLUMN 6 WILL ALWAYS BE 
ZEROES TO ALLOW FOR THE CALCULATION OF OVER 85 
VALUES. AC ACCUMULATES THE AIR CONDITIONER 
DATA. 

ALLGR : REAL ARRAY USED TO STORE GROWTH RATES OF THE 

OCTANTS OF THE CONTOURS 

AVGAGE : A 3-DIMENSIONAL ARRAY WHOSE ROWS CORRESPOND TO 

THE 6 AGE CATEGORIES; COLUMNS CORRESPOND TO THE 
8 OCTANTS AND PLANES CORRESPOND TO THE LDN 
VALUES. 

AVGFIN : A 3-DIMENSIONAL ARRAY WHOSE ROWS CORRESPOND TO 

THE 7 FAMILY INCOME CATEGORIES; COLUMNS 
CORRESPOND TO THE 8 OCTANTS AND PLANES 
CORRESPOND TO THE LDN VALUES. 

AVGHMV : A 3-DIMENSIONAL ARRAY WHOSE ROWS CORRESPOND TO 

THE 7 HOME VALUE CATEGORIES; COLUMNS CORRESPOND 
TO THE 8 OCTANTS AND PLANES CORRESPOND TO THE 
LDN VALUES. 

BANDGRO : REAL ARRAY REPRESENTING BAND GROWTH RATES WITH 

THE FINAL 8 OCTANT VARIABLES REPRESENTING THE 
INNER CONTOUR (8 OCTANTS, 5 BANDS, WHERE LAST 
BAND IS "OVER XX" 

DIF : TEMPORARY VARIABLE USED TO ACCUMULATE 

DIFFERENCES THAT ARE MULTIPLIED BY A FACTOR. 

DIFHHS : TEMPORY USED FOR DIFFERENCING THE NO. OF 

HOUSEHOLDS. 
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c 
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DIF1 : 
DIF2 : 
DIF3 : 
DIF4 : 

DIF5 : 
DIF6 : 
DIF7 : 
DIFS : 
CMS IT Y : 
DWL1 : 

EDI 6 : 


FACT : 


FACTHM : 
FACTIN : 
FPOP : 
GRATE : 

GRIOCT : 
GROTH : 
GR9 : 


HHS : 

HMVALT : 

I : 

IALLPOP: 

ICT : 

IDIF : 


TEMPORARY FOR DIFFERENCING THE AVG. HOME VALUES 
TEMPORARY FOR DIFFERENCING THE YRS EDUCATION. 
TEMPORARY FOR DIFFERENCING THE OCCUPATIONS. 

TEMPORARY FOR DIFFERENCING THE YRS EDUCATION 

OVER 16. 

TEMPORARY FOR DIFFERENCING THE MANAGER/PROF. 
TEMPORARY FOR DIFFERENCING THE SINGLE DWELLINGS 
TEMPORARY FOR DIFFERENCING THE AC UNITS. 

TEMPORARY FOR DIFFERENCING THE TVS. 

TEMPORARY FOR COMPUTING POPULATION DENSITY. 
ACCUMULATION ARRAY FOR SINGLE DWELLING DATA; 

ROWS ARE OCTANTS, COLUMNS ARE LDN VALUES. 
3-DIMENSIONAL ARRAY WHICH ACCUMULATES THE 
EDUCATIONAL CATEGORIES; ROW 1 IS 16YRS+ AND 
ROW 2 IS ALL CATEGORIES, COLUMNS ARE 
OCTANTS, PLANES ARE LDN VALUES. 

ARRAY OF MULTIPLYING FACTORS USED FOR COMPUTING 
THE AVERAGE AGE DATA; FOR EXAMPLE, THE 1ST 
ELEMENT IS 19 WHICH CORRESPONDS TO THE AVERAGE 
AGE OF TOE 1ST AGE CATEGORY (18-20) BEING 
SUMMED FOR AVERAGING. 

SIMILAR TO FACT ABOVE BUT USED FOR AVERAGING 
THE HOME VALUES. 

SIMILAR TO FACT ABOVE BUT USED FOR AVERAGING 
THE FAMILY INCOMES. 

TEMPORARY USED FOR DIFFERENCING THE FAMILY 
POPULATION. 

REAL VARIABLE REPRESENTING THE GROWTH RATE, 

EITHER EXTRACTED FROM ARRAY BANDGRO, OR 
USER-DESIGNATED; THIS VARIABLE IS THE GROWTH 
RATE USED TO PROJECT A CORRESPONDING POPULATION 
REAL VARIABLE REPRESENTING THE GROWTH RATE FOR 
A PARTICULAR OCTANT 

REAL ARRAY USED TO STORE THE GROWTH RATES FOR 

AN OCTANT WITH MULTIPLE POLYGONS 

REAL VARIABLE, REPRESENTING THE GROWTH RATE FOR 

THE "ALL OCTANTS" CASE, COMPUTED USING THE ORIGINAL 

CUMULATIVE POPULATION AND THE GROWN CUMULATIVE 

POPULATION 

TEMPORARY USED FOR ACCUMULATING THE HOUSEHOLD 
DATA; ROWS ARE OCTANTS, COLUMNS ARE LDN VALUES. 
TEMPORARY USE FOR TOTALING THE HOME VALUES. 

GENERAL PURPOSE DO LOOP INDEX. 

INTEGER ARRAY USED TO STORE POPULATION VALUES OF 
THE OCTANTS OF THE CONTOURS, NEEDED TO COMPUTE 
BAND GROWTH RATES 

COUNTER USED AS INDEX TO ACCUMULATE GROWTH RATES 
AND AND POPULATIONS IN SINGLE OCTANT 
INTEGER VARIABLE USED FOR TEMPORARY STORAGE TO 
REPRESENT THE NUMBER OF PEOPLE FOR A PARTICULAR 
OCTANT AND A PARTICULAR CONTOUR BAND 
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IGRPOP 

IOCT2 

IOCTSUM 

ITEMP 

J 

JJ 

JJJ 

JPTR 

K 

KK 

KDUM1 

KDUM2 

LEW 

LOC 

N 

POP 

POPDEN 


POPIMP : 

POPO : 

POP1 : 

PRATIO : 
SAREA : 
SHHS : 
SHMV 
SPOP : 
SPOPO : 


SPOP70 : 
STIN : 


INTEGER ARRAY USED TO STORE THE POPULATIONS OF 
AN OCTANT WITH MULTIPLE POLYGONS 
INTEGER VARIABLE USED TO IDENTIFY OCTANT, THUS 
CAN ACCOUNT FOR MULTIPLE POLYGONS IN 1 OCTANT 
INTEGER VARIABLE REPRESENTING THE TOTAL VALUE 
OF THE POPULATION IN A GIVEN OCTANT 
INTEGER VARIABLE USED TO INDEX INTO ARRAY IALLPOP 
TO FIND POPULATION OF INNER CONTOUR 
DO LOOP INDEX FOR 600 LOOP( OCTANT INDEX ) . 

INDEX USED INTO ARRAY ALINE, USED TO STORE 
AREA INFORMATION 

INDEX INTO ARRAY ALINE, USED TO STORE. AREA 

INFORMATION 

LDN VALUE INDEX. 

GENERAL PURPOSE DO LOOP INDEX 

INDEX INTO ARRAY POPDEN, USED TO EXTRACT AREA 

INFORMATION 

DUMMY PARAMETER FOR GETCON ROUTINE. 

DUMMY PARAMETER FOR GETCON ROUTINE. 

INTEGER REPRESENTATION FOR THE CURRENT LDN 
VALUE. 

INTEGER VARIABLE USED TO LOCATE STORAGE OF 

CONTOUR AREAS, ON AN OCTANT BY OCTANT BASIS 

INTEGER VARIABLE REPRESENTING THE NUMBER 

OF YEARS OF POPULATION UPDATE 

TEMPORARY USED FOR DIFFERENCING THE POPULATION 

DENSITIES. 

3-DIMENSIONAL ARRAY USED FOR ACCUMULATING THE 
1981 POPULATIONS AND CORRESPONDING AREAS; COL 1 
1981 POPULATION, COL 2 IS AREA, ROWS ARE 
OCTANTS AND COLUMNS ARE LDNS. 

ARRAY FOR ACCUMULATING THE IMPACTED POPULATIONS 
ROWS ARE OCTANTS AND COLUMNS ARE LDNS. 

REAL VARIABLE USED FOR DIFFERENCING POPULATION 
VALUES 

REAL VARIABLE USED TO STORE 'THE POPULATION OF A 
POLYGON WITHIN AN OCTANT OF A CONTOUR 
TEMPORARY FOR CALCULATING THE POPULATION RATIOS 
TEMPORARY FOR SIM OF AREAS. 

TEMPORARY FOR SUM OF HOUSEHOLDS. 

TEMPORARY FOR SUM OF HOME VALUES. 

TEMPORARY USED FOR SUMMING POPULATIONS. 

REAL VARIABLE USED TO STORE SUMMED POPULATIONS; 
CONTAINS THE SUMMED POPULATIONS WITHOUT APPLI- 
CATION OF THE GROWTH RATE 

TEMPORARY USED FOR SUMMING 1970 POPULATIONS. 
TEMPORARY USED FOR SUMMING TOTAL INCOME. 
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SWTSUM : TEMPORARY USED FOR SUMMING WEIGHTED SUMS FOR 

THE AVERAGE AGES. 

S18T65 : SUM OF TOE 18 TO 65 AGE GROUPS 

TERM : REAL VARIABLE USED TO COMPUTE THE WEIGHTED 

GROWTH RATE OF AN OCTANT IN THE CASE OF AN 
OCTANT WITH MULTIPLE POLYGONS 

TFMPOP : ARRAY FOR ACCUMULATING FAMILY POPULATIONS; ROWS 

ARE OCTANTS, COLUMNS ARE LDNS. 

TINCOM : TEMPORARY FOR SUMMING THE TOTAL INCOMES. 

TMNPRO : 3-DIMENSIONAL ARRAY FOR ACCUMULATING THE 

OCCUPATIONS; ROW 1 IS MANAGER/PROF ONLY AND 
ROW 2 IS ALL CATEGORIES; COLS ARE OCTANTS 
PLANES ARE LDNS. 

TPOP70 : TEMPORARY FOR ACCUMULATING THE 1970 POPULATIONS 

ROWS ARE OCTANTS, COLUMNS ARE LDNS. 

TV : TEMPORARY FOR ACCUMULATING TV DATA. ROWS ARE 

OCTANTS, COLUMNS ARE LDNS. 

T10 ; TEMPORARY FOR SUMMING OCCUPATIONAL CATEGORIES. 

T18T65 : TEMPORARY FOR TOTALING THE AVERAGE AGE DATA. 

T9 : TEMPORARY FOR TOTALING THE EDUCATIONAL DATA. 

WT : TEMPORARY FOR TOTALING THE WEIGHTED SIMS OF THE 

FAMILY INCOME. 

WTSUM : GENERAL PURPOSE TEMPORARY FOR STORAGE OF 

WEIGHTED SUMS. 

WTSUM1 : SAME AS WTSUM. 

WT1 : TEMPORARY FOR SUMMING THE WTSUMl WEIGHTED SUMS. 

X : TEMPORARY VARIABLE USED IN LOGARITHMIC 

CALCULATIONS. 


COMMON VARIABLES: 

' /COMM/ 

SEE SUBROUTINE RDREPT COMMON BLOCK DESCRIPTION. 


A- 5 7 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE HEADER 


OBJECTIVE: 

TO PRINT OCTAtfT REPORT HEADERS 


PARAMETERS: 

INPUT: 

KOCT : THE OCTANT NUMBER OF THE OCTANT BEING 

PRINTED. ( 1 - 8 ) 

NAME : THE AIRPORT NAME FOR THE HEADING ( 80 R1 ) 

KADRS : THE AIRPORT ADDRESS ( 80 R1 ) 

. OUTPUT: 

ICASE : 60-WORD ARRAY CONTAINING CASE NAME 

MINUS : INTEGER VARIABLE, NON-ZERO FOR LANDUSE 

INSULATION COUNTERMEASURES, THAT REPRE- 
SENTS THE DB DECREMENT FOR THE 65 LDN 
(AND GREATER) VALUE 

CALLING ROUTINE: 

AUREPT 


SUBPROGRAMS USED: 

GETCON : COMPASS SUBROUTINE THAT SAMPLES REGISTER VALUES 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 


IOCT 

: 3 LETTER MNEUMONIC FOR OCTANT 

AS 

FOLLOWS: 


1 = ESE 

2 = SSE 3 = SSW 

4 

= WSW 5 = 


6 = NNW 

7 = NNE 8 = ENE 



IR1 

: INTEGER 

VALUE OF REGISTER 1 



IR2 

: INTEGER 

VALUE OF REGISTER 2 



IR3 

: INTEGER 

VALUE OF REGISTER 3 



K 

: GENERAL 

PURPOSE IMPLIED DO LOOP 

INDEX. 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

TAPE 6 - LINE PRINTER FILE WITH CARRIAGE CONTROL 
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SUBROUTINE CHART 
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OBJECTIVE: 

TO PRINT THE 16 DEMOGRAPHIC VARIABLES BY LDN BAND VALUES 
IN REPORT FORM USING THE LINE PRINTER. 

PARAMETERS: 

INPUT: 

TABLE : THE 3-DIMENSIONAL TABLE STORING THE 

DEMOGRAPHIC VARIABLES TO BE PRINTED. 

LDNS : A VECTOR STORING THE LDN VALUES IN 

ASCENDING ORDER. 

N : THE OCTANT INDEX TO TABLE 

NCON : STARTING INDEX TO THE LDN TABLE. IF NCON 

IS GREATER THAN ONE THEN 1 OR MORE LDN 
CONTOURS WERE MISSING. 

IYEAR : USER-DESIGNATED YEAR OF POPULATION UPDATE 

GR2 : ' USER-DESIGNATED GROWTH RATE 


CALLING ROUTINE: 
AUREPT 


SUBPROGRAMS USED: 

BREAK : BREAK BREAKS OUT 5 COLUMNS OUT OF THE 

3-DIMENSIONAL TABLE TO A VECTOR CALLED LINE. 
CHART1 : ALTERNATIVE OCTANT PRINTING ROUTINE USED WHEN 

LDN CONTOURS ARE MISSING (NCON GT 1) 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

I : IMPLIED DO LOOP INDEX. 

J : IMPLIED DO LOOP INDEX. 

K : DO LOOP INDEX FOR REPLACING VERY SMALL AREA 

QUANTITIES WITH A . 

LINE : GENERAL PURPOSE 5 WORD VECTOR USED TO ENCODE 

BAND VALUES FOR CONVENIENT PRINTING. 

COMMON VARIABLES: 

NONE 

FILES USED: 

OUTPUT: 

TAPE6 : SEQUENTIAL LINE PRINTER FILE FOR OCTANT 

REPORT. 
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SUBROUTINE CHART1 


C 
C 
C 
C 

C OBJECTIVE: 

C TO PRINT OCTANT REPORTS TO THE PRINTER AS AN ALTERNATIVE 

C TO SUBROUTINE CHART WHEN LDN CONTOUR ARE MISSING. 

C 

C 


c 

PARAMETERS: 




c 

INPUT: 




c 

TABLE 

• 

• 

3-DIMENS IONAL TABLE STORING THE 


c 



DEMOGRAPHIC VARIABLES. 


c 

LDNS 

• 

6 WORD VECTOR STORING THE RANGE OF 

LDN 

c 



CONTOUR VALUES. 


c 

N 

• 

THE. OCTANT NUMBER TO BE PRINTED. 


c 

NCON 

: 

THE STARTING INDEX IN THE TABLE OF 

LDNS 

c 



WHICH COMPENSATES FOR MISSING CONTOURS. 

c 

I YEAR 

: 

USER-DESIGNATED YEAR OF POPULATION 

UPDATES 

c 

GR2 

; 

USER-DESIGNATED GROWTH RATE 


c 

OUTPUT: 




c 

c 

NONE 




c 

c 

CALLING ROUTINE: 




c 

c 

CHART 




c 

c 

SUBPROGRAMS USED: 





C BREAK : BREAK BREAKS OUT 5 COLUMNS OUT OF THE 

C 3-DIMENSIONAL TABLE TO A VECTOR CALLED LINE. 

C 
C 

C FUNCTIONS USED: 

C NONE 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LOCAL VARIABLE DICTIONARY: 


I 

INDEX 


INT 

J 

JCON 

K 

KCON 

LINE 


IMPLIED DO LOOP INDEX. 

ARRAY INDEX TO LDN TABLE USED TO COMPUTE 
MISSING CONTOUR VALUES AT LOW END OF CONTOUR 
RANGE. 

THE INTERVAL BETWEEN CONTOUR VALUES; EXPECTED 
VALUE IS 5. 

IMPLIED DO LOOP INDEX. 

THE NUMBER OF CONTOUR BANDS THAT WERE FOUND. 

DO LOOP INDEX USED TO CHECK FOR SMALL AREAS. 

THE NUMBER OF CONTOUR BANDS THAT WERE NOT FOUND 
GENERAL PURPOSE 5 WORD VECTOR USED TO ENCODE 
BAND VALUES FOR CONVENIENT PRINTING. 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 

OUTPUT: 

SEQUENTIAL LINE PRINTER FILE FOR OCTANT REPORT. 
SCRATCH: 

NONE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE BREAK 


OBJECTIVE: 

BREAK CREATES A VECTOR FROM A SELECTED ROW OF THE 
3-DIMENSIONAL TABLE. 


PARAMETERS: 

INPUT:. 

TABLE : THE 3-DIMENSIONAL TABLE OF DEMOGRAPHIC 

VARIABLES. 

KROW : THE SELECTED ROW. TO BE VECTORIZED. 

KPLANE : THE SELECTED PLANE WHERE THE SELECTED ROW 

RESIDES. 

OUTPUT: 

LINE : THE VECTOR REMOVED FROM THE 3-DIMENSIONAL 

TABLE. 


CALLING ROUTINE: 

CHART, CHART1 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

I : DO LOOP INDEX USED FOR COLUMNS OF TABLE. 


COMMON VARIABLES: 
NONE 


FILES USED: 

NONE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE SUMARY 


OBJECTIVE: 

SUMARY COMPUTES AND PRINTS A NOISE IMPACT SUMMARY REPORT 
BASED ON THE POPULATION DATA AND THE LEN VALUES. 


METHOD: 

SIX VALUES ARE CALCULATED FOR EACH OCTANT AND A TOTAL 

FOR EACH VALUE IS COMPUTED FOR ALL OCTANTS. THE VALUES 

ARE NAMED AND COMPUTED AS FOLLOWS: 

IMPACTED POPULATION THE SUM OF THE 1981 POPULATIONS 

FOR ALL BANDS OF AN OCTANT. 

LEVEL 'WEIGHTED POP. THE POPULATION EXPOSED TO A 

CERTAIN LDN BAND IS WEIGHTED. 

THE FOLLOWING FORMULA IS USED TO 
COMPUTE THE WEIGHT CORRESPONDING 
TO A GIVEN LDN VALUE: 

.103LDN 
.000124(10 ) 

•03LDN .08LDN 

.2(10 )+. 000143(10 ) 

THE LEVEL WEIGHTED POPULATION IS 
THE SUM OF THE PRODUCTS OF THE 
POPULATIONS AND THEIR WEIGHTS. 

NOISE IMPACT INDEX THE LEVEL WEIGHTED POPULATION 

DIVIDED BY THE TOTAL POPULATION. 

EQUIVALENT NOISE LEVEL A TABLE OF WEIGHTS IS CALCULATED 

FOR EACH LDN VALUE USING THE 
FORMULA GIVEN ABOVE. THE LDN 
VALUES ARE THE BAND AVERAGES; EG 
65 - 70 HAS VALUE 67. THE ENL IS 
OBTAINED BY INTERPOLATING THE 
Nil IN THE TABLE OF WEIGHTS. 

ONCE A WEIGHT IS KNOWN A 
CORRESPONDING LDN CAN BE 
INTERPOLATED. THIS LDN IS THE 
ENL. 

PERCENT HIGHLY ANNOYED THE PRODUCT OF THE Nil AND THE 

CONSTANT 36.9. 

POP. HIGHLY ANNOYED THE PRODUCT OF THE LWP AND THE 

CONSTANT .369 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PARAMETERS: 

INPUT: 

ALINE : 2-DIMENSIONAL ARRAY CONTAINING AREAS 

OF 65, 75, 85 DB CONTOURS. 

ICASE : 60-WORD ARRAY CONTAINING CASE TITLE 

MINUS : INTEGER VALUE, NON-ZERO FOR LANDUSE 

INSULATION COUNTERMEASURES, REPRE- 
SENTING , FOR THOSE CASES, THE DB 
DECREMENT FOR LDN VALUES GREATER 
THAN OR EQUAL TO 65 LDN 

TABLE : THE 3-DIMENSIONAL TABLE USED HERE TO 

SUPPLY THE POPULATION DATA. 

LDNS : 8 WORD VECTOR CONTAINING THE CONTOUR LDN 

VALUES. 

LINE1 : 80 WORD ARRAY FOR STORAGE OF THE AIRPORT 

NAME IN 80 R1 FORMAT. 

LINE2 : 80 WORD ARRAY FOR STORAGE OF THE AIRPORT 

ADDRESS IN 80 R1 FORMAT. 

OUTPUT: 

NONE 


CALLING ROUTINE: 
AUREPT 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

INTERP : USED TO FIND A WEIGTING FACTOR GIVEN AN Nil; 

THEN INTERPOLATES AN LDN FROM THE WEIGHT. 

PERCNT : USED TO CALCULATE THE WEIGHTS USING THE FORMULA 

DESCRIBED ABOVE. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LOCAL VARIABLE DICTIONARY : 

DBVAL : INTEGER VARIABLE USED TO REPRESENT CONTOUR 

VALUE WHEN PRINTING OUT AREA OF CONTOUR 

ENL : A 9 WORD VECTOR USED TO STORE THE COMPUTED 

EQUIVALENT NOISE LEVEL FOR PRINTING BY OCTANT 
AND TOTAL OF ALL OCTANTS. 

I : GENERAL PURPOSE DO LOOP INDEX. 

IDASH : WORD USED TO STORE THE SYMBOL TO 

BE USED AS SUBSTITUTE VALUE FOR Nil, ENL, 

PCTHA, AND POPHA WHEN A CORRESPONDING 
BASELINE CASE POPULATION VALUE EQUALS ZERO 
(UTILIZED ONLY IN COUNTERMEASURE CASES) . 

IR1 : VALUE OF REGISTER 1 

IR2 : VALUE OF REGISTER 2 

IR3 : VALUE OF REGISTER 3 

J : DO LOOP INDEX FOR OCTANTS. 

JJ : INDEX USED FOR DECODING Nil, ENL, PCTHA, 

AND POPHA 

K : INDIRECT DO LOOP INDEX BASED ON J TO REVERSE 

INDEX THE ARRAYS. 

LDNTAB : A 5 WORD VECTOR USED TO STORE THE LDN BAND 

AVERAGE VALUES; 67,72,77,82,85 
THE LAST VALUE IS FOR OVER 85 AREA. 

LINE : 9-WORD ARRAY USED IN PRINTING OF Nil, ENL, 

PCTHA, AND POPHA 

LWP : 9 WORD VECTOR USED TO STORE THE COMPUTED LEVEL 

WEIGHTED POPULATIONS BY OCTANT AND TOTAL. 

Nil : 9 WORD VECTOR USED TO STORE THE COMPUTED NOISE 

IMPACT INDEX BY OCTANT AND TOTAL. 

PCTHA : 9 WORD VECTOR USED TO STORE THE PERCENT HIGHLY 

ANNOYED BY OCTANT AND TOTAL. 

POPHA : 9 WORD VECTOR USED TO STORE THE POPULATION 

HIGHLY ANNOYED BY OCTANT AND TOTAL. 

TOTPOP : USED TO SUM THE RUNNING TOTAL BY OCTANTS. 

WTS : A 5 WORD VECTOR USED TO STORE THE WEIGHTING 

FACTORS FOR THE LDN BANDS. 

XLDN : USED FOR FLOATING PT. STORAGE OF THE INTEGER 

VALUES IN THE ARRAY LDNTAB. 


COMMON VARIABLES: 
NONE 


FILES USED: 

INPUT: 

NONE 


C OUTPUT: 

C TAPE6 - SEQUENTIAL - LINE PRINTER FILE FOR SUMMARY 

C REPORT. 
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FUNCTION INTERP 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OBJECTIVE: 

TO RETURN AN EQUIVALENT NOISE LEVEL VALUE IN LDN FORM 
GIVEN A NOISE IMPACT INDEX. 


METHOD: 

INTERP USES THE GIVEN NOISE IMPACT INDEX TO FIND AN ENTRY 
IN THE WEIGHTS TABLE THAT IS GREATER THAN THE NOISE 
IMPACT INDEX. THE WEIGHT TABLE STORES THE WEIGHTS IN 
ASCENDING ORDER CORRESPONDING TO THE LDN BAND AVERAGE 
VALUES. INTERP THEN COMPUTES THE PROPORTION OF TOE NOISE 
IMPACT INDEX BETWEEN THE WEIGHT ENTRIES IN THE TABLE. 

THIS PROPORTION IS THEN USED TO COMPUTE THE APPROPRIATE 
LDN WHICH IS PROPORTIONATELY LOCATED BETWEEN THE TWO 
CORRESPONDING LDN ENTRIES IN THE LDN TABLE. IF THE NOISE 
IMPACT INDEX IS LESS THAN THE FIRST ENTRY IN THE WEIGHT 
TABLE, THE FLOOR VALUES FOR WEIGHTS AND LDNS ARE ZEROES. 
IF NO ENTRY IN THE WEIGHTS TABLE IS GREATER THAN THE 
NOISE IMPACT INDEX THEN INTERP RETURNS THE GREATEST LDN 
VALUE IN THE LDN TABLE AS THE EQUIVALENT NOISE LEVEL. 


PARAMETERS: 

INPUT: 

LDNS : A 6 WORD VECTOR STORING THE CONTOUR LDN 

VALUES IN ASCENDING ORDER. 

WTS : A 6 WORD VECTOR OF WEIGHTING FACTORS 

CORRESPONDING TO THE ELEMENTS OF THE 
VECTOR LDNS. 

Nil : THE GIVEN NOISE IMPACT INDEX FOR WHICH AN 

EQUIVALENT NOISE LEVEL IS TO BE RETURNED. 

OUTPUT: 

NONE 


CALLING ROUTINE: 
SUMARY 


SUBPROGRAMS USED: 
NONE 


FUNCTIONS USED: 

FLOAT : INLINE FUNCTION FOR INTEGER TO REAL CONVERSION. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LOCAL VARIABLE DICTIONARY: 

BOT : THE FLOOR VALUE OF THE TWO WEIGHTS VALUES THAT 

ENCLOSE THE NOISE IMPACT INDEX. 

I : DO LOOP INDEX FOR WEIGHTS AND LDNS TABLES. 

INTERP : THE FUNCTION VARIABLE USED TO RETURN THE ENL. 

XBOT : THE FLOOR VALUE OF THE TWO LDN VALUES THAT 

ENCLOSE THE ENL. 

XLDNS : THE CEILING VALUE OF THE TWO LDN VALUES THAT 

ENCLOSE THE ENL. 
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SUBROUTINE LASTPAG 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OBJECTIVE: 

TO WRITE A SEPARATE REPORT CONSISTING OF TOE 
NOISE IMPACT SUMMARY REPORT 


PARAMETERS: 

NONE 


CALLING ROUTINE: 
AUREPT 


SUBPROGRAMS USED: 
NONE 


LOCAL VARIABLE DICTIONARY: 

IPAGE : INTEGER ARRAY USED TO READ IN AND WRITE OUT REPORT 
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FUNCTION PERCNT 


C 

c 
c 
c 
c 

C OBJECTIVE: 

C TO CADCULATE A WEIGHTING FACTOR GIVEN AN LDN VALUE. 

C 
C 

C METHOD: 

C 
C 
C 

Q 

c 
c 
c 

C .103LDN 

C PERCENT HIGHLY ANNOYED = .000124(10 ) 

C 

C . 0 3LDN .08LDN 

C .2(10 )+. 000143(10 ) 

C 

C 

C PARAMETERS: 

C INPUT: 

C LDN : THE LON VALUE USED TO CADCULATE A WEIGHT. 

C OUTPUT: 

C NONE 

C 

C 

C CALLING ROUTINE: 

C SUMARY 

C 

C 

SUBPROGRAMS USED: 

NONE 


PERCNT IMPLEMENTS THE FOLLOWING FORMULA WHICH IS THE 
CURVE OF PERCENT OF POPULATION "HIGHLY ANNOYED" AS A 
FUNCTION OF NOISE LEVEL IN LDN. THE PERCENT HIGHLY 
ANNOYED CURVE IS NORMALIZED TO UNITY AT 75LDN BY 
DIVIDING THE PERCENT HIGHLY ANNOYED BY 36.9, THE 
PERCENT OF POPULATION HIGHLY ANNOYED AT 75LDN. 


C 


c 

FUNCTIONS USED: 


c 

NONE 


v.. 

C 

0 

U1CAL VARIABLE DIC 

T IONARY : 

C 

El : 

TEMPORARY FOR VALUE OF .103*LDN 

0 

E2 : 

TEMPORARY FOR VALUE OF .03*LDN 

c 

E3 : 

TEMPORARY FOR VALUE OF .08*LDN 

c 

PERCNT : 

THE FUNCTION VARIABLE USED TO RETURN THE 

c 


COMPUTED WEIGHT. 
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APPENDIX. B 


SITE II DATA BASE INSTALLATION 
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Four tapes were delivered by CACI, Inc. containing the 1981 SITE II 
database. These tapes contain data reflecting the 1970 census with 1980 
to 1981 updates. The population data was contained on tapes NASA13, 
NASA14 and NASA15. The zip code data was contained on tape NASA00. 
NASA13 contained 15 PRU's or blocks of data each block 4424 characters 
in length, NASA14 and NASA15 contained 18 blocks each 4424 characters in 
length. The setup programs SETUP 15 and SETUP18 were used to read the 
tape files according to the above mentioned specifications and write an 
intermediate file TAPE7 to be processed by the conversion prog ran COMMN. 
The data from these tapes were converted from IBM character format to 
CDC CYBER binary data and written to the NOS labelled tapes ND0775 and 
ND1062. 

Figures 16, 17 and 18 contain the job stream of control statements 
used to convert the IBM formatted population data and generate the NOS 
tapes. Figures 19, 20 and 21 illustrate the source for SETUP15, SETUP18 
and CONMN, respectively. Figure 22 illustrates the job stream used to 
read the tape NASA00, convert the data using progran CONZP and generate 
the new labelled tape ND1226. Figure 23 contains the source code for 
program CONZP. 
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CONMN, T300 f 0160000. 

USER. 

CHARGE. 

FTN, I=SETUP15. 

FILE (TAPE5, RT=F, BT=K, FL=4424, MBL=4424, RB=1) 

VSN, TAPE5=NASA13. 

REQUEST, TAPE5, F=S, NT, D=1600, CV=EB, PO=R, LB=KU . N0NLRC=NASA1 3 
LDSET, FILES=TAPE5. 

LGO. 

RETURN, TAPE5. 

REWIND, TAPE7, LGO. 

GET, CONMN. 

FTN, I=CONMN. 

LOAD (LGO) 

NOGO(CONMS) 

VSN, TAPE1=ND0775. 

LABEL, TAPE1, FI=AL, SI=MN1, QN=1, W, PO=W. 

CONMS. 

REWIND, TAPES . 

COPYBF, TAPES, TAPE1 . 

RETURN, TAPES. 

LABEL, TAPE1, FI=AK, SI=MN1, QN=9999 . 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE 1. 

RETURN, TAPE3. 
o 
o 
o 

LABEL, TAPE1, FI=IN, SI=MN1, QN=9999. 

CONMS. 

REWIND, TAPES. 

COPYBF, TAPES, TAPE 1. 

RETURN, TAPES. 


Figure 16. -Step 1 job for the SITE II database conversion. 
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C0NM2, T300, CM60000. 

USER. 

CHARGE. 

GET, SETUP18. 

FTN, I=SETUP18. 

FILE(TAPE5, RT=F, BT=K, FL=4424, MBL=4424, RB=1) 

VSN, TAPE5=NASA14 . 

REQUEST ( TAPE5 , F=S, NT, E>=1GOO, LB=KU, CV=EB, POR) N0NLRC=NASA14 
LDSET ( FILES=TAPE5) 

DGO. 

RETURN, TAPE5. 

REWIND, TAPE7,LG0. 

GET, CONMN. 

FTN, I=CONMN . 

LOAD (DGO) 

NOGO (CONMS) 

VSN(TAPE1=ND0775) 

LABEL(TAPE1, FI=IA, SI=MN1, QN=9999, PO=W) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8. 

LAB EL ( TAPE 1, FI=KS, SI=MN1, QN=9999) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAP El . 

RETURN, TAPE8 . 
o 
o 

LABEL(TAPE1, FI=MS, SI=MN1, QN=9999) 

CONMS. 

REWIND, TAPE8 . 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8 . 

RETURN, TAPE1. 

VSN, TAPE1=ND1062 . 

LABEL(TAPE1, FI=MO, SI=MN2, QN=1, W, POW) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8 . 
o 
o 

LABEL (TAPE1, FI=NY, SI=MN2, QN=9999) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8 . 

Figure 17. -Step 2 job for the SITE II database conversion. 
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C0NM3, T300, CM60000 . 

USER. 

CHARGE. 

GET, SETUP 18 . 

FILE(TAPE5, RT=F, BT=K, FL=4424, MBL=4424, RB=1) 

VSN (TAPE5=NASA15) 

REQUEST (TAPE5, F=S, NT, D=1600, LB=KU, CV=EB, POR) N0NLRC=NASA15 
LDS ET ( F I LES=TAPE 5 ) 

LGO. 

RETURN, TAPE5. 

REWIND, TAPE7, LGO. 

GET, CONMN. 

FTN, I=CONMN . 

LOAD (LGO) 

NOGO(CONMS) 

VSN, TAPE1=ND1062 . 

LABEL(TAPE1, FI=NC, SI=MN2, QN=9999) 

CONMS . 

REWIND, TAPES. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8. 

LABEL (TAPE1, FI=ND, SI=MN2, QN=9999) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8. 
o 
o 

LABEL(TAPE1, FI=WY, SI=MN2, QN=9999) 

CONMS. 

REWIND, TAPE8. 

COPYBF, TAPE8, TAPE1 . 

RETURN, TAPE8 . 


Figure 18. -Step 3 job for the SITE II database conversion 
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PROGRAM SETUP1 5 (OUTPUT, TAPE5=/4424, TAPE7=/4424) 
DIMENSION IN (443) 

IREC=1 

INREC=0 

IFICT=1 

INP=0 

IOT=0 

5 READ(5, 100) IN 

IF(EOF(5) ) 11, 7 
7 INP=INEM-1 

WRITE (7, 100) IN 
IOT>=IOT+l 

IF(INP.GE.2) INREOINREC+6 
GO TO 5 

11 ENDFILE 7 

IF(IFICT.EQ.15)GO TO 10 

PRIWT 130, IREC, I OT, INREC, INP 

IREC=1 

INREC=0 

IFICT=IFICT+1 

INP=0 

GO . TO 5 

10 PRINr 120, IOT 

STOP, 

100 FORMAT (150A10, 150A10, 142A10, A4) 

120 FORMAT (33H1NORMAL END OF PROCESSING, COUNT=,I5) 

110 FORMAT(8X,9X, 19) 

130 FORMAT (25H0ENDFILE WRITTEN. COUNTS=, 415) 

END 


Figure 19. -Setup program for 15 records. 
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PROGRAM SETUP 18 (OUTPUT, TAPE5=/4424, TAPE7=/4424) 
DIMENSION IN (443) 

IREC=1 

INREC=0 

IFICT=1 

INP=0 

IOT=0 

5 READ (5, 100) IN 

IF(EOF(5) ) 11, 7 
7 INP=INP+1 

WRITE (7, 100) IN 
IOT=IOT+l 

IF(INP.GE.2) INREC=INREC+6 
GO TO 5 

11 ENDFILE 7 

IF(IFICT.EQ.18)GO TO 10 

PRINr 130, IREC, IOT, INREC, INP 

IREC=1 

INREC=0 

IFICT=IFICT+1 

INP=0 

GO TO 5 

10 PRIOT 120, IOT 

STOP 

100 FORMAT (150A10, 150A10, 142A10, A4) 

120 FORMAT ( 3 3H1NORMAL END OF PROCESSING, COUNT=, 15) 

110 FORMAT (3X,9X, 19) 

130 FORMAT ( 2 5H0ENDF I LE WRITTEN. COUNTS=,4I5) 

ENO 


Figure 20. -Setup program to process 18 records. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


5 


7 


11 . 


10 

100 

110 

120 


PROGRAM CONMN (OUTPUT, TAPE6=OUTPUT, TAPE7=/4424, TAPE8) 

PROGRAM TO CONVERT MAIN DATABASE FROM CHAR TO BINARY NASA/LRC . 
DATABASE CREATED ON IBM 370 WITH FORMAT 

BLKSIZE=4424, LRECL=4424, DEN=3, STANDARD EBCDIC 
PARAMS 

INPUT REC = 736 CHARACTERS 

INPUT BLOCK = 4416 CHARACTERS 

INPUT BLOCK FACTOR = 6 

OUTPUT REC = 89 BINARY WORDS 3 DBL PRECISION=92E 

OUTPUT BLOCK = 534 BINARY WORDS 3 DBL PRECISION=552E 

OUTPUT BLOCK FACTOR = 6 

DIMENSION INDEX (800) , 11(552) , FOREC1 (92, 6) , FOREC2(46, 6) 

COMMON IOREC1 (92, 6) 

DOUBLE PRECISION FOREC2 
EQUIVALENCE (II (1) , IOREC1 (1, 1) ) 

EQUIVALENCE ( IOREC1 (1,1), FOREC1 (1,1)) 

EQUIVALENCE ( IOREC1 ( 1, 1) , FOREC2 ( 1, 1) ) 

IREC=1 

INREC=0 

IFILE=7 

INP=0 

IOT=l 

CALL OPENMS(8, INDEX, 800, 0) 

READ ( IFILE, 100) ( ( IOREC1 ( I, J) , 1=1, 2) , (FOREC1 ( I, J) , 1=3, 4 ) , 

1 (FOREC2(I, J) , 1=3, 5) , (I0REC1(I, J) , 1=11, 92) , J=l, 6) 

IF(EOF( IFILE)) 10,7 
INP=INP+-1 

IF(INP.GE.2) IREC=I0REC1 (2, 1) 

CALL WRnWS(8, 11,552, IOT) 

IOT=IOT+l 

IF( INP.GE.2) INREC=INREC+6 

GO TO 5 

IOT^lOT-1 

PRINT 120, INP, IOT, INREC, IREC 
CALL CLOSMS (8) 

STOP 

PRINT 110, INP 
STOP 

FORMAT(8X, 6(219, 2F7.4, 3F16.0, 8218) ) 

FORMAT ( 3 3HONORMAL END OF PROCESSING, COUNT=,I5) 

FORMAT (140FILE WRITTEN., 12H INPUT BLOCKS, 15, 13HOUTPOT BLOCKS, 

1, 15, 13H INPUT RECORDS, 15, 7HCOUNTER, 15) 

END 


Figure 21. -Database conversion program. 
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CONZ IP, T300, CM60000 . 

USER. 

CHARGE. 

GET, CONZ P. 

FTN, I=CONZP. 

FILE(TAPE7, RT=F, BT=K, FL=3208, MBL=3208, RB=1) 
VSN, TAPE7=NASA00 . 

REQUEST (TAPE7, F=S, NT, 0=1600, LB=KU, CV=EB, PO=R) 
LD5 ET ( F I LES=TAPE7 ) 

LGO. 

RETURN, TAPE7. 

VSN, TAPE2=ND1225. 

LABEL (TAPE2, SI=ZIPCDS, FI=ZIPALL, QN=1, W, PO=W) 
COPYBF, TAPES , TAPE 2 . 

RETURN, TAPE 2. 


Figure 22. -Job to process the zip code portion of SITE II. 
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PROGRAM CONZP (OUTPUT, TAPE6=OUTPUT r TAPE7=/3208, TAPE8) 


C 

C PROGRAM TO CONVERT ZIP DATABASE FROM CHARACTER TO 
C BINARY FOR NASA/LRC. 

C DATABASE CREATED ON IBM 370 WITH FORMAT: 

C LRECL=3204/ BLKSIZE=3208, RECFM=VS, DEN=3 

C PARAMS: 

C INPUT RECORD = 32 CHARACTERS 

C INPUT BLOCK = 3200 CHARACTERS 

C OUTPUT BLOCK FACTOR = 100 

C OUTPUT RECORD = 5 BINARY WORDS 

C OUTPUT BLOCK =500 BINARY WORDS 

C OUTPUT BLOCK FACTOR = 100 

C 

DIMENSION INDEX(4000) , 11(500) 

DIMENSION IOREC ( 5, 100) , FOREC(5, 100) 

EQUIVALENCE (1 1(1) , IOREC (1, 1) ) 

EQUIVALENCE ( IOREC ( 1, 1) , FOREC (1,1)) 

I0T=1 

INP=0 

CALL OPENMS(8, INDEX, 4000, 0) 

5 READ(7, 100) ( (IOREC(I, J) , 1=1, 3) , FOREC(4, J) , IOREC(5, J) , J=l, 100) 
IF(EOF(7) ) 10,7 
7 INP=INPfl 

CALL WRITMS(8, II, 500, IOT) 

IOT=IOT+l 
GO TO 5 

10 PRINT 110, INP 
STOP 

100 FORMAT(8X, 100 (15, 219, F5.4, 14) ) 

110 FORMAT (33HNORMAL END OF PROCESSING, COUNT=,I5) 

END 


Figure 23. -Conversion program to process zip code data 
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